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INTRODUCTION TO DIFFERENTIAL 


EQUATIONS 


1.1 Definitions and Terminology 

1.2  Initial- Value Problems 

1.3 Differential Equations as Mathematical Models 
CHAPTER 1 IN REVIEW 


The words differential and equations certainly suggest solving some kind of 
equation that contains derivatives y’, y”,.... Analogous to a course in algebra and 
trigonometry, in which a good amount of time is spent solving equations such as 
xX + 5x + 4 = 0 for the unknown number x, in this course one of our tasks will be 
to solve differential equations such as y” + 2y’ + y = 0 for an unknown function 
y = $x). 

The preceding paragraph tells something, but not the complete story, about the 
course you are about to begin. As the course unfolds, you will see that there is more 
to the study of differential equations than just mastering methods that someone has 
devised to solve them. 

But first things first. In order to read, study, and be conversant in a specialized 
subject, you have to learn the terminology of that discipline. This is the thrust of the 
first two sections of this chapter. In the last section we briefly examine the link 
between differential equations and the real world. Practical questions such as How 
fast does a disease spread? How fast does a population change? involve rates of 
change or derivatives. As so the mathematical description—or mathematical 


model—of experiments, observations, or theories may be a differential equation. 
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CHAPTER 1 


INTRODUCTION TO DIFFERENTIAL EQUATIONS 


1.1 


DEFINITIONS AND TERMINOLOGY 


REVIEW MATERIAL 


Definition of the derivative 

Rules of differentiation 

Derivative as a rate of change 

First derivative and increasing/decreasing 
Second derivative and concavity 


INTRODUCTION The derivative dy/dx of a function y = #(x) is itself another function d'G 


found by an appropriate rule. The function y = e 
by the Chain Rule its derivative is dy/dx = 0.2xe°'“” If we replace e 


0.1x° ig differentiable on the interval (—%, ©), and 


0.1" on the right-hand side of 


the last equation by the symbol y, the derivative becomes 


dy 
— = 0.2xy. 1 
ES xy (1) 


Now imagine that a friend of yours simply hands you equation (1)—you have no idea how it was 
constructed —and asks, What is the function represented by the symbol y? You are now face to face 
with one of the basic problems in this course: 


How do you solve such an equation for the unknown function y = (x)? 


A DEFINITION The equation that we made up in (1) is called a differential 
equation. Before proceeding any further, let us consider a more precise definition of 
this concept. 


| DEFINITION 1.1.1 Differential Equation 


An equation containing the derivatives of one or more dependent variables, 
with respect to one or more independent variables, is said to be a differential 
equation (DE). 


To talk about them, we shall classify differential equations by type, order, and 
linearity. 


CLASSIFICATION BY TYPE Ifan equation contains only ordinary derivatives of 
one or more dependent variables with respect to a single independent variable it is 
said to be an ordinary differential equation (ODE). For example, 


A DE can contain more 
than one dependent variable 


dy dy dy dx dy 
SE + 6y = 0, di AFS 
a E di a” ii Ar dt 


=2x+y (2) 


are ordinary differential equations. An equation involving partial derivatives of 
one or more dependent variables of two or more independent variables is called a 
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partial differential equation (PDE). For example, 


u Pu u gu ðu ðu ðv 
EE SE and eg (3) 
dr du ox dr dr oy Ox 


are partial differential equations.” 

Throughout this text ordinary derivatives will be written by using either the 
Leibniz notation dy/dx, d*y/dx*, d*y/dx*, . . . or the prime notation y’, y”, y”,.... 
By using the latter notation, the first two differential equations in (2) can be written 
a little more compactly as y’ + 5y = e* and y” — y’ + 6y = 0. Actually, the prime 
notation is used to denote only the first three derivatives; the fourth derivative is 
written y instead of vi". In general, the nth derivative of y is written d"y/dx" or y™®. 
Although less convenient to write and to typeset, the Leibniz notation has an advan- 
tage over the prime notation in that it clearly displays both the dependent and 
independent variables. For example, in the equation 


unknown function 
or dependent variable 


ag det 
dt? 


E independent variable 


it is immediately seen that the symbol x now represents a dependent variable, 
whereas the independent variable is t£. You should also be aware that in physical 
sciences and engineering, Newton’s dot notation (derogatively referred to by some 
as the “flyspeck” notation) is sometimes used to denote derivatives with respect 
to time 1 Thus the differential equation d*s/dt? = —32 becomes ï = —32. Partial 
derivatives are often denoted by a subscript notation indicating the indepen- 
dent variables. For example, with the subscript notation the second equation in 
(3) becomes Uy, = Uy — Ju 


CLASSIFICATION BY ORDER The order of a differential equation (either 
ODE or PDE) is the order of the highest derivative in the equation. For example, 


second order sl T first order 
d?y dy\3 F 
ES + (| — Au =e 


is a second-order ordinary differential equation. First-order ordinary differential 
equations are occasionally written in differential form M(x, y) dx + N(x, y) dy = 0. 
For example, if we assume that y denotes the dependent variable in 
(y — x) dx + 4xdy = 0, then y’ = dy/dx, so by dividing by the differential dx, we 
get the alternative form 4xy’ + y = x. See the Remarks at the end of this section. 

In symbols we can express an nth-order ordinary differential equation in one 
dependent variable by the general form 


F(x, y, y',.--,y¥”) = 0, (4) 


where F is a real-valued function of n + 2 variables: x, y, y’, . . . , vil For both prac- 
tical and theoretical reasons we shall also make the assumption hereafter that it is 
possible to solve an ordinary differential equation in the form (4) uniquely for the 


“Except for this introductory section, only ordinary differential equations are considered in A First 
Course in Differential Equations with Modeling Applications, Ninth Edition. In that text the 

word equation and the abbreviation DE refer only to ODEs. Partial differential equations or PDEs 
are considered in the expanded volume Differential Equations with Boundary- Value Problems, 
Seventh Edition. 
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INTRODUCTION TO DIFFERENTIAL EQUATIONS 


highest derivative y™ in terms of the remaining n + 1 variables. The differential 
equation 


n 


J. ' (n—1) 5 
ES: Jy wn Area YTO (5) 
where fis a real-valued continuous function, is referred to as the normal form of (4). 
Thus when it suits our purposes, we shall use the normal forms 


dy dy 

— =f, and —=fay.y’ 

ES fœ, y) ES fF, y, y’) 

to represent general first- and second-order ordinary differential equations. For example, 
the normal form of the first-order equation 4xy’ + y = xis y’ = (x — y)/4x; the normal 
form of the second-order equation y” — y’ + 6y = Ois y” = vi — 6y. See the Remarks. 


CLASSIFICATION BY LINEARITY An nth-order ordinary differential equation (4) 
is said to be linear if F is linear in y, y’,... , y. This means that an nth-order ODE is 
linear when (4) is an(x)y + ayy"? + +++ + ay + ao(x)y — g(x) = 0 or 


n Dal 


a(x) E T a,- 1%) > 
dx" 


d 
Sech eee SS + ag(x)y = g(x). (6) 


Two important special cases of (6) are linear first-order (n = 1) and linear second- 
order (n = 2) DEs: 
dy *y dy 
oul + ao(x)y = gx) = and a(x) (5 + A(X) Tt awy = g(x). O) 
dx dx dx 


In the additive combination on the left-hand side of equation (6) we see that the char- 
acteristic two properties of a linear ODE are as follows: 


e The dependent variable y and all its derivatives y’, y",..., y are of the 
first degree, that is, the power of each term involving y is 1. 

e The coefficients ao, a), ..., an of y, y’,..., y™® depend at most on the 
independent variable x. 


The equations 


(y — x)dx + 4xdy=0, y"—2y’+y=0, and +x 5y = ei 


are, in turn, linear first-, second-, and third-order ordinary differential equations. We 
have just demonstrated that the first equation is linear in the variable y by writing it in 
the alternative form 4xy’ + y = x. A nonlinear ordinary differential equation is sim- 
ply one that is not linear. Nonlinear functions of the dependent variable or its deriva- 
tives, such as sin y or ei, cannot appear in a linear equation. Therefore 


nonlinear term: nonlinear term: nonlinear term: 
coefficient depends on y nonlinear function of y power not 1 
| m | any | 
, X : y 
1—y)y'+2y=e, —5 tsiny=0, and =, +y’=0 
(l= yp y ES y gi? 


are examples of nonlinear first-, second-, and fourth-order ordinary differential equa- 
tions, respectively. 


SOLUTIONS As was stated before, one of the goals in this course is to solve, or 
find solutions of, differential equations. In the next definition we consider the con- 
cept of a solution of an ordinary differential equation. 
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DEFINITION 1.1.2 Solution of an ODE 


Any function œ, defined on an interval J and possessing at least n derivatives 
that are continuous on 7, which when substituted into an nth-order ordinary 
differential equation reduces the equation to an identity, is said to be a 
solution of the equation on the interval. 


In other words, a solution of an nth-order ordinary differential equation (4) is a func- 
tion dh that possesses at least n derivatives and for which 


F(x, (x), dit), ..., dito) = 0 for all x in I. 


We say that dh satisfies the differential equation on /. For our purposes we shall also 
assume that a solution ¢ is a real-valued function. In our introductory discussion we 
saw that y = e®! is a solution of dy/dx = 0.2xy on the interval (—%, ©). 

Occasionally, it will be convenient to denote a solution by the alternative 
symbol y(x). 


INTERVAL OF DEFINITION You cannot think solution of an ordinary differential 
equation without simultaneously thinking interval. The interval J in Definition 1.1.2 
is variously called the interval of definition, the interval of existence, the interval 
of validity, or the domain of the solution and can be an open interval (a, b), a closed 
interval [a, b], an infinite interval (a, ©), and so on. 


l EXAMPLE 1 Verification of a Solution 


Verify that the indicated function is a solution of the given differential equation on 
the interval (—%, ©), 


(a) dy/dx = xy'?; y=} (b) y" — 2y' + y=0; y=xe* 


SOLUTION One way of verifying that the given function is a solution is to see, after 
substituting, whether each side of the equation is the same for every x in the interval. 


(a) From 


: 1 
left-hand side: —=—(4-x) =-x, 


1 1/2 1 1 
right-hand side: xy” =x. Le SA ( e) =—x, 


we see that each side of the equation is the same for every real number x. Note 
that y"? = ix? is, by definition, the nonnegative square root of 4% x°. 


(b) From the derivatives y’ = xe* + e* and y” = xe” + 2e* we have, for every real 


number x, 
left-hand side: y” — 2y' + y = (xe + 2e”) — 2(xe* + e*) + xæ = 0, 
right-hand side: 0. a 


Note, too, that in Example 1 each differential equation possesses the constant so- 
lution y = 0, -~ < x < œ. A solution of a differential equation that is identically 
zero on an interval / is said to be a trivial solution. 


SOLUTION CURVE The graph of a solution dh of an ODE is called a solution 
curve. Since ¢ is a differentiable function, it is continuous on its interval J of defini- 
tion. Thus there may be a difference between the graph of the function and the 
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T T i T e 
(a) function y = 1/x, x #0 
ya 
1 = D 
+— i H S 


(b) solution y = 1/x, (0, œ) 


FIGURE 1.1.1 


The function y = 1/x 
is not the same as the solution y = 1/x 


graph of the solution dh. Put another way, the domain of the function dh need not be 
the same as the interval 7 of definition (or domain) of the solution dh. Example 2 
illustrates the difference. 


| EXAMPLE 2 Function versus Solution 


The domain of y = 1/x, considered simply as a function, is the set of all real num- 
bers x except 0. When we graph y = 1/x, we plot points in the xy-plane corre- 
sponding to a judicious sampling of numbers taken from its domain. The rational 
function y = 1/x is discontinuous at 0, and its graph, in a neighborhood of the ori- 
gin, is given in Figure 1.1.1(a). The function y = 1/x is not differentiable at x = 0, 
since the y-axis (whose equation is x = 0) is a vertical asymptote of the graph. 
Now y = 1/x is also a solution of the linear first-order differential equation 
xy’ + y = 0. (Verify.) But when we say that y = 1/x is a solution of this DE, we 
mean that it is a function defined on an interval J on which it is differentiable and 
satisfies the equation. In other words, y = 1/x is a solution of the DE on any inter- 
val that does not contain 0, such as (—3, —1), (4, 10), (—œ, 0), or (0, ©). Because 
the solution curves defined by y = 1/x for —3 < x < —1 and 5 <x < 10 are sim- 
ply segments, or pieces, of the solution curves defined by y = 1/x for ~% <x <0 
and 0 < x < ™, respectively, it makes sense to take the interval J to be as large as 
possible. Thus we take / to be either (—™, 0) or (0, ©). The solution curve on (0, œ) 
is shown in Figure 1.1.1(b). D 


EXPLICIT AND IMPLICIT SOLUTIONS You should be familiar with the terms 
explicit functions and implicit functions from your study of calculus. A solution in 
which the dependent variable is expressed solely in terms of the independent 
variable and constants is said to be an explicit solution. For our purposes, let us 
think of an explicit solution as an explicit formula y = d(x) that we can manipulate, 
evaluate, and differentiate using the standard rules. We have just seen in the last two 
examples that y = xx, y =xe*, and y= 1/x are, in turn, explicit solutions 
of dy/dx = xy", vi — 2y' + y = 0, and xy’ + y = 0. Moreover, the trivial solu- 
tion y = 0 is an explicit solution of all three equations. When we get down to 
the business of actually solving some ordinary differential equations, you will 
see that methods of solution do not always lead directly to an explicit solution 
y = (x). This is particularly true when we attempt to solve nonlinear first-order 
differential equations. Often we have to be content with a relation or expression 


G(x, y) = 0 that defines a solution dh implicitly. 


| DEFINITION 1.1.3 Implicit Solution of an ODE 


A relation G(x, y) = 0 is said to be an implicit solution of an ordinary 
differential equation (4) on an interval 7, provided that there exists at least 
one function ¢ that satisfies the relation as well as the differential equation 
on J. 


It is beyond the scope of this course to investigate the conditions under which a 
relation G(x, y) = 0 defines a differentiable function dh. So we shall assume that if 
the formal implementation of a method of solution leads to a relation G(x, y) = 0, 
then there exists at least one function œ that satisfies both the relation (that is, 
G(x, ġ(x)) = 0) and the differential equation on an interval 7. If the implicit solution 
G(x, y) = 0 is fairly simple, we may be able to solve for y in terms of x and obtain 
one or more explicit solutions. See the Remarks. 
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Y4 
2 | EXAMPLE 3 Verification of an Implicit Solution 
z The relation x? + y? = 25 is an implicit solution of the differential equation 
x dy x 
Sa 8 
~ (8) 
on the open interval (—5, 5). By implicit differentiation we obtain 
(a) implicit solution d d d d 
SW EEGENEN ET 
x +y =25 dx dx dx dx 


Solving the last equation for the symbol dy/dx gives (8). Moreover, solving 
x? + y? =25 for y in terms of x yields y = +\/25 — x. The two functions 
y = ġ(x) = V25 — x and y = h(x) = — V25 — x satisfy the relation (that is, 
x? + bj = 25 and x? + 3 = 25) and are explicit solutions defined on the interval 
(—5, 5). The solution curves given in Figures 1.1.2(b) and 1.1.2(c) are segments of 
the graph of the implicit solution in Figure 1.1.2(a). a 


Any relation of the form x? + y? — e = 0 formally satisfies (8) for any constant c. 
However, it is understood that the relation should always make sense in the real number 


(b) explicit solution system; thus, for example, if c = —25, we cannot say that x? + y? + 25 =0 is an 
EE implicit solution of the equation. (Why not?) 
yı 7 =R aa X 


Because the distinction between an explicit solution and an implicit solution 
should be intuitively clear, we will not belabor the issue by always saying, “Here is 
an explicit (implicit) solution.” 


FAMILIES OF SOLUTIONS The study of differential equations is similar to that of 
integral calculus. In some texts a solution dh is sometimes referred to as an integral 
of the equation, and its graph is called an integral curve. When evaluating an anti- 
derivative or indefinite integral in calculus, we use a single constant c of integration. 
Analogously, when solving a first-order differential equation F(x, y, y’) = we 
usually obtain a solution containing a single arbitrary constant or parameter c. A 


(c) explicit solution solution containing an arbitrary constant represents a set G(x, y, c) = 0 of solutions 
called a one-parameter family of solutions. When solving an nth-order differential 

yy = —-V25—-x7, 5 <x <5 : P ' Gë 8 ` . 
equation F(x, y, y',...,y™)=0, we seek an n-parameter family of solutions 
FIGURE 1.1.2 An implicit solution G(x, Y, C1, Co... Cn) = 0. This means that a single differential equation can possess 
and two explicit solutions of y' = —x/y an infinite number of solutions corresponding to the unlimited number of choices 


for the parameter(s). A solution of a differential equation that is free of arbitrary 
parameters is called a particular solution. For example, the one-parameter family 
y = cx — x cos x is an explicit solution of the linear first-order equation xy’ — y = 
x? sin x on the interval (—%, oc). (Verify.) Figure 1.1.3, obtained by using graphing soft- 
ware, shows the graphs of some of the solutions in this family. The solution y = 
—x cos x, the blue curve in the figure, is a particular solution corresponding to c = 0. 
Similarly, on the interval (—%, ©), y = cje” + c2xe* is a two-parameter family of solu- 
tions of the linear second-order equation y” — 2y' + y = 0 in Example 1. (Verify.) 
Some particular solutions of the equation are the trivial solution y = 0 (c1 = c2 = 0), 
y = xe* (cy = 0, c2 = 1), y = Se* — 2xe* (c1 = 5, c2 = —2), and so on. 

Sometimes a differential equation possesses a solution that is not a member of a 
family of solutions of the equation—that is, a solution that cannot be obtained by spe- 
F IG URE 1 : 1 3 Some solutions of cializing any of the parameters in the family of solutions. Such an extra solution is called 
ae a ST a singular solution. For example, we have seen that y = Je 7 and y = O are solutions of 

the differential equation dy/dx = xy"? on (—%, ~). In Section 2.2 we shall demonstrate, 
by actually solving it, that the differential equation dy/dx = xy"? possesses the one- 
parameter family of solutions y = (I r+ DH When c = 0, the resulting particular 
solution is y = 5 x*. But notice that the trivial solution y = 0 is a singular solution, since 
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Ya 


c=-l, 
x<0 


(b) piecewise-defined solution 


FIGURE 1.1.4 Some solutions of 
xy’ — 4y =0 


it is not a member of the family y = (i Ke cht: there is no way of assigning a value to 
the constant c to obtain y = 0. 

In all the preceding examples we used x and y to denote the independent and 
dependent variables, respectively. But you should become accustomed to seeing 
and working with other symbols to denote these variables. For example, we could 
denote the independent variable by ¢ and the dependent variable by x. 


l EXAMPLE 4 Using Different Symbols 


The functions x = cı cos 4t and x = c3 sin 4t, where cı and c3 are arbitrary constants 
or parameters, are both solutions of the linear differential equation 


x" + 16x = 0. 
For x = cı cos 4t the first two derivatives with respect to ¢ are x’ = —4c; sin 4t 
and x” = — 16c; cos 4t. Substituting x” and x then gives 


x" + 16x = —16c, cos 4t + 16(c, cos 41) = 0. 
In like manner, for x = cz sin 4t we have x” = —16c2 sin 4t, and so 
x" + 16x = —16c,sin 4t + 16(c, sin An = 0. 


Finally, it is straightforward to verify that the linear combination of solutions, or the 
two-parameter family x = cı cos 4t + c2 sin 4t, is also a solution of the differential 
equation. EI 


The next example shows that a solution of a differential equation can be a 
piecewise-defined function. 


l EXAMPLE ZS A Piecewise-Defined Solution 


You should verify that the one-parameter family y = cx* is a one-parameter family 


of solutions of the differential equation xy’ — 4y = 0 on the inverval (—™, zl, See 
Figure 1.1.4(a). The piecewise-defined differentiable function 


=x, x<0 
e ae oa 
A. x20 


is a particular solution of the equation but cannot be obtained from the family 
y = cxf by a single choice of c; the solution is constructed from the family by choos- 
ing c = —1 for x <Oandc = 1 for x = See Figure 1.1.4(b). E 


SYSTEMS OF DIFFERENTIAL EQUATIONS Up to this point we have been 
discussing single differential equations containing one unknown function. But 
often in theory, as well as in many applications, we must deal with systems of 
differential equations. A system of ordinary differential equations is two or more 
equations involving the derivatives of two or more unknown functions of a single 
independent variable. For example, if x and y denote dependent variables and t 
denotes the independent variable, then a system of two first-order differential 
equations is given by 


d 

e = TL, 3 

i (9) 
= = g(t, x, y). 
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A solution of a system such as (9) is a pair of differentiable functions x = dutt, 
y = datt, defined on a common interval J, that satisfy each equation of the system 
on this interval. 


l REMARKS 


(i) A few last words about implicit solutions of differential equations are in 
order. In Example 3 we were able to solve the relation x? + y? = 25 for 
y in terms of x to get two explicit solutions, h(x) = V25 — x? and 
(x) = — V25 — x’, of the differential equation (8). But don’t read too much 
into this one example. Unless it is easy or important or you are instructed to, 
there is usually no need to try to solve an implicit solution G(x, y) = 0 for y 
explicitly in terms of x. Also do not misinterpret the second sentence following 
Definition 1.1.3. An implicit solution G(x, y) = 0 can define a perfectly good 
differentiable function ¢ that is a solution of a DE, yet we might not be able to 
solve G(x, y) = 0 using analytical methods such as algebra. The solution curve 
of dh may be a segment or piece of the graph of G(x, y) = 0. See Problems 45 
and 46 in Exercises 1.1. Also, read the discussion following Example 4 in 
Section 2.2. 


(ii) Although the concept of a solution has been emphasized in this section, 
you should also be aware that a DE does not necessarily have to possess 
a solution. See Problem 39 in Exercises 1.1. The question of whether a 
solution exists will be touched on in the next section. 


(iii) It might not be apparent whether a first-order ODE written in differential 
form M(x, y)dx + N(x, y)dy = 0 is linear or nonlinear because there is nothing 
in this form that tells us which symbol denotes the dependent variable. See 
Problems 9 and 10 in Exercises 1.1. 


(iv) It might not seem like a big deal to assume that F(x, y, y’,..., vy) = 0 can 
be solved for y®, but one should be a little bit careful here. There are exceptions, 
and there certainly are some problems connected with this assumption. See 
Problems 52 and 53 in Exercises 1.1. 


(v) You may run across the term closed form solutions in DE texts or in 
lectures in courses in differential equations. Translated, this phrase usually 
refers to explicit solutions that are expressible in terms of elementary (or 
familiar) functions: finite combinations of integer powers of x, roots, exponen- 
tial and logarithmic functions, and trigonometric and inverse trigonometric 


functions. 

(vi) If every solution of an nth-order ODE F(x, y, y', . . . y) = 0 on an inter- 
val J can be obtained from an n-parameter family G(x, y, C1, C2, ..., Cn) = 0 by 
appropriate choices of the parameters c; i = 1, 2, ..., n, we then say that the 


family is the general solution of the DE. In solving linear ODEs, we shall im- 
pose relatively simple restrictions on the coefficients of the equation; with these 
restrictions one can be assured that not only does a solution exist on an interval 
but also that a family of solutions yields all possible solutions. Nonlinear ODEs, 
with the exception of some first-order equations, are usually difficult or impos- 
sible to solve in terms of elementary functions. Furthermore, if we happen to 
obtain a family of solutions for a nonlinear equation, it is not obvious whether 
this family contains all solutions. On a practical level, then, the designation 
“general solution” is applied only to linear ODEs. Don’t be concerned about 
this concept at this point, but store the words “general solution” in the back of 
your mind—we will come back to this notion in Section 2.3 and again in 
Chapter 4. 
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INTRODUCTION TO DIFFERENTIAL EQUATIONS 


EXERCISES 1.1 


Answers to selected odd-numbered problems begin on page ANS-1. 


In Problems 1-8 state the order of the given ordinary differ- 
ential equation. Determine whether the equation is linear or 
nonlinear by matching it with (6). 


1. (1 — x)y” — 4xy'+ 5y = cos x 


a dy\* 
x us ( d +y=0 
dx? dx 


p 


3: py — Py" +6y=0 


d? d 
4. E EE EE 
dr? dr 
d? 2 
sell (2) 
dx? dx 
PE o E 
` dt? R? 


7. (sin 0)y” — (cos 0)y' = 2 
-2 
dit Cl +x=0 


In Problems 9 and 10 determine whether the given first-order 
differential equation is linear in the indicated dependent 
variable by matching it with the first differential equation 
given in (7). 


9. (y? — 1) dx + x dy = 0; in y; in x 
10. u dv + (v + uv — ue") du = 0; in v; in u 
In Problems 11-14 verify that the indicated function is an 


explicit solution of the given differential equation. Assume 
an appropriate interval J of definition for each solution. 


11. 2y +y=0; y=e*? 


dy 6 
12. — + 20y = 24; == 
dt E S 5 


13. y” — 6y' + 13y = 0; y= e** cos 2x 


14. y” + y= tanx; y= —(cos x)In(sec x + tan x) 


In Problems 15-18 verify that the indicated function 
y = (x) is an explicit solution of the given first-order 
differential equation. Proceed as in Example 2, by consider- 
ing dh simply as a function, give its domain. Then by consid- 
ering ¢ as a solution of the differential equation, give at least 
one interval / of definition. 


15. (y - y =y-—x+8 y=x+4Vx+2 


16. y! =25 +y% y= 5 tan 5x 

17. y! =2xy*s y=1/(4- x’) 

18. 2y' = y? cosx; y= (1 — sin x)! 

In Problems 19 and 20 verify that the indicated expression is 
an implicit solution of the given first-order differential equa- 
tion. Find at least one explicit solution y = ġ(x) in each case. 


Use a graphing utility to obtain the graph of an explicit solu- 
tion. Give an interval / of definition of each solution dh, 


SEH 


19 = (Xx — Id — 2X); 1 ( 
k = ; In 
t A 


d 


20. 2xy dx + (x? —y)dy=0; -2x y +y =1 

In Problems 21 24 verify that the indicated family of func- 
tions is a solution of the given differential equation. Assume 
an appropriate interval J of definition for each solution. 


dP ? 
21. Ën P); p=—* 
dt 1+ ce 


d oo ae : 
22. L 2xy = 1; y= f e dt +cje* 
dx 0 


x 


23. + 4y =0; y= ce” + cxe 
x 


d d? d 
24. GE 12x; 
x x x 


y= ax! + ox + coxlnx + 4° 


25. Verify that the piecewise-defined function 


is a solution of the differential equation xy’ — 2y = 0 
on (—%, %). 


26. In Example 3 we saw that y = d(x) = V25 — x° and 
y = bx) = —V25 — x are solutions of dy/dx = 
—x/y on the interval (—5, 5). Explain why the piecewise- 
defined function 


8 V25 — x7, äere 
yO |-V25—- 2, 05x <5 


is not a solution of the differential equation on the 
interval (—5, 5). 
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In Problems 27-30 find values of m so that the function 


IX + 


y = e™ is a solution of the given differential equation. 


27. vi + 2y=0 
29. y” — 5y' + by =0 


28. Sy’ = 2y 
30. 2y” + Ty’ — 4y =0 


In Problems 31 and 32 find values of m so that the function 
y = x” is a solution of the given differential equation. 


31. xy" + 2y'=0 
32. xy" — Ixy’ + 15y =0 


In Problems 33—36 use the concept that y = c, ~% < x < ©, 
is a constant function if and only if y’ =0 to determine 
whether the given differential equation possesses constant 
solutions. 


33. 3xy' + 5y = 10 

34. y! =y?+2y-3 

35. y- 1)y'=1 

36. y” + 4y’ + 6y = 10 

In Problems 37 and 38 verify that the indicated pair of 


functions is a solution of the given system of differential 
equations on the interval (—%, 2%). 


dx d'r 
7. —=x+ . — = 4y + e 
3 e x + 3y 38 Te yte 
dy d’y 
be" bdt Se GZ A Se 


x= eo + 3e%, x = cos 2t + sin 2t + Zei, 


cos 2t — sin 2t b 


e 


y=-e 7 + Set y= i 


Discussion Problems 


39. Make up a differential equation that does not possess 
any real solutions. 


40. Make up a differential equation that you feel confident 
possesses only the trivial solution y = 0. Explain your 
reasoning. 


41. What function do you know from calculus is such that 
its first derivative is itself? Its first derivative is a 
constant multiple k of itself? Write each answer in 
the form of a first-order differential equation with a 
solution. 


42. What function (or functions) do you know from calcu- 
lus is such that its second derivative is itself? Its second 
derivative is the negative of itself? Write each answer in 
the form of a second-order differential equation with a 
solution. 


1.1 DEFINITIONS AND TERMINOLOGY ° 11 


43. Given that y = sin x is an explicit solution of the 


d 
first-order differential equation SS = V1 — y’. Find 
x 


an interval J of definition. [Hint: I is not the interval 


Een 


44. Discuss why it makes intuitive sense to presume that 
the linear differential equation y” + 2y’ + 4y = 5 sin t 
has a solution of the form y = A sin t + B cos t, where 
A and B are constants. Then find specific constants A 
and B so that y = A sin t + B cos t is a particular solu- 
tion of the DE. 


In Problems 45 and 46 the given figure represents the graph 
of an implicit solution G(x, y) = 0 of a differential equation 
dy/dx = f(x, y). In each case the relation G(x, y) = 0 
implicitly defines several solutions of the DE. Carefully 
reproduce each figure on a piece of paper. Use different 
colored pencils to mark off segments, or pieces, on each 
graph that correspond to graphs of solutions. Keep in mind 
that a solution ¢ must be a function and differentiable. Use 
the solution curve to estimate an interval J of definition of 
each solution dh, 


45. ya 


FIGURE 1.1.5 Graph for Problem 45 


46. ya 


FIGURE 1.1.6 Graph for Problem 46 


47. The graphs of members of the one-parameter family 
x? + y? = 3cxy are called folia of Descartes. Verify 
that this family is an implicit solution of the first-order 
differential equation 


dy _ y? — 2x) 
dx xy — x) 
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48. 


49. 


50. 


51. 


52. 


53. 


54. 


55. 


° CHAPTER 1 


The graph in Figure 1.1.6 is the member of the family of 
folia in Problem 47 corresponding to c = 1. Discuss: 
How can the DE in Problem 47 help in finding points 
on the graph of x? + y? = 3xy where the tangent line 
is vertical? How does knowing where a tangent line is 
vertical help in determining an interval J of definition 
of a solution dh of the DE? Carry out your ideas, 
and compare with your estimates of the intervals in 
Problem 46. 


In Example 3 the largest interval J over which the 
explicit solutions y = ġı(x) and y = @o(x) are defined 
is the open interval (—5, 5). Why can’t the interval Z of 
definition be the closed interval [—5, 5]? 


In Problem 21 a one-parameter family of solutions of 
the DE P’ = P(1 — P) is given. Does any solution 
curve pass through the point (0, 3)? Through the 
point (0, 1)? 


Discuss, and illustrate with examples, how to solve 
differential equations of the forms dy/dx = f(x) and 
Ste 


The differential equation x(y')* — 4y’ — 12x = 0 has 
the form given in (4). Determine whether the equation 
can be put into the normal form dy/dx = f(x, y). 


The normal form (5) of an nth-order differential equa- 
tion is equivalent to (4) whenever both forms have 
exactly the same solutions. Make up a first-order differ- 
ential equation for which F(x, y, y’) = 0 is not equiva- 
lent to the normal form dy/dx = f(x, y). 


Find a linear second-order differential equation 
F(x, y, y', y”) = 0 for which y = cx + c2x? is a two- 
parameter family of solutions. Make sure that your equa- 
tion is free of the arbitrary parameters cı and c2. 


Qualitative information about a solution y = (x) of a 
differential equation can often be obtained from the 
equation itself. Before working Problems 55-58, recall 
the geometric significance of the derivatives dy/dx 
and d?y/dx?. 


x 


Consider the differential equation dy/dx ze". 
(a) Explain why a solution of the DE must be an 
increasing function on any interval of the x-axis. 

(b) What are lim- dy/dx and lim dy/dx? What does 

this peges about a sditon curve as x > £0? 


(c) Determine an interval over which a solution curve is 
concave down and an interval over which the curve 
is concave up. 

(d) Sketch the graph of a solution y = (x) of the dif- 
ferential equation whose shape is suggested by 
parts (a)—(c). 


56. 


57. 


58. 


INTRODUCTION TO DIFFERENTIAL EQUATIONS 


Consider the differential equation dy/dx = 5 — y. 


(a) Either by inspection or by the method suggested in 
Problems 33-36, find a constant solution of the DE. 


(b) Using only the differential equation, find intervals on 
the y-axis on which a nonconstant solution y = d(x) 
is increasing. Find intervals on the y-axis on which 
y = (x) is decreasing. 


Consider the differential equation dy/dx = y(a — by), 
where a and b are positive constants. 


(a) Either by inspection or by the method suggested 
in Problems 33—36, find two constant solutions of 
the DE. 


(b) Using only the differential equation, find intervals on 
the y-axis on which a nonconstant solution y = (x) 
is increasing. Find intervals on which y = (x) is 
decreasing. 


(c) Using only the differential equation, explain why 
y = a/2b is the y-coordinate of a point of inflection 
of the graph of a nonconstant solution y = (x). 


(d) On the same coordinate axes, sketch the graphs of 
the two constant solutions found in part (a). These 
constant solutions partition the xy-plane into three 
regions. In each region, sketch the graph of a non- 
constant solution y = dc) whose shape is sug- 
gested by the results in parts (b) and (c). 


Consider the differential equation y’ = y? + 4. 


(a) Explain why there exist no constant solutions of 
the DE. 

(b) Describe the graph of a solution y = d(x). For 
example, can a solution curve have any relative 
extrema? 


(c) Explain why y = 0 is the y-coordinate of a point of 
inflection of a solution curve. 

(d) Sketch the graph of a solution y = d(x) of the 
differential equation whose shape is suggested by 
parts (a)—-(c). 


Computer Lab Assignments 


In Problems 59 and 60 use a CAS to compute all derivatives 
and to carry out the simplifications needed to verify that the 
indicated function is a particular solution of the given differ- 
ential equation. 


59. 


60. 


y® — 20y" + 158y” — 580y’ + 841y = 0; 
y = xe™ cos 2x 
xy" + Zivil + 20xy' — 78y = 0; 

— om cos(5 In x) 3 sin(5 In x) 


xX 
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solutions of the DE 


> 
x 


FIGURE 1.2.1 ` Solution of 
first-order IVP 


solutions of the DE 


m=y, 
L > 
= = ZS 


FIGURE 1.2.2 Solution of 
second-order IVP 


where yo, Yi... 


y' (xo) = yh... 
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1.2 INITIAL-VALUE PROBLEMS 

REVIEW MATERIAL 
e Normal form of a DE 
e Solution of a DE 
e Family of solutions 
INTRODUCTION We are often interested in problems in which we seek a solution y(x) of a 
differential equation so that y(x) satisfies prescribed side conditions—that is, conditions imposed on 
the unknown y(x) or its derivatives. On some interval J containing xo the problem 

d"y 

= seca e 
du (1) 
Subject to: y@o) = Yo Yo) = Yr eV" PCO) = Yn- 


»Yn-1 are arbitrarily specified real constants, is called an initial-value 
problem (IVP). The values of y(x) and its first n — 1 derivatives at a single point xo, vol = yo, 
, y" D(x9) = yn—1, are called initial conditions. 


FIRST- AND SECOND-ORDER IVPS The problem given in (1) is also called an 
nth-order initial-value problem. For example, 


dy 


Solve: = = E E 
dx (2) 
Subject to: y(%) = Yo 
dy 
d SS E ; 
an Solve E fy, y') (3) 


Subject to: Mt = Yo. Co) = Yı 


are first- and second-order initial-value problems, respectively. These two problems 
are easy to interpret in geometric terms. For (2) we are seeking a solution y(x) of the 
differential equation y’ = f(x, y) on an interval 7 containing xo so that its graph passes 
through the specified point (xo, yo). A solution curve is shown in blue in Figure 1.2.1. 
For (3) we want to find a solution y(x) of the differential equation y” = f(x, y, y’) on 
an interval J containing xo so that its graph not only passes through (xo, yo) but the slope 
of the curve at this point is the number yı. A solution curve is shown in blue in 
Figure 1.2.2. The words initial conditions derive from physical systems where the 
independent variable is time £ and where y(fo) = yo and y'(fo) = vu represent the posi- 
tion and velocity, respectively, of an object at some beginning, or initial, time fo. 

Solving an nth-order initial-value problem such as (1) frequently entails first 
finding an n-parameter family of solutions of the given differential equation and then 
using the n initial conditions at x9 to determine numerical values of the n constants in 
the family. The resulting particular solution is defined on some interval J containing 
the initial point xo. 


l EXAMPLE 1 Two First-Order IVPs 


In Problem 41 in Exercises 1.1 you were asked to deduce that y = ce* is a one- 
parameter family of solutions of the simple first-order equation y’ = y. All the 
solutions in this family are defined on the interval (—~, ~). If we impose an initial 
condition, say, y(0) = 3, then substituting x = 0, y = 3 in the family determines the 
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FIGURE 1.2.3 


Solutions of two IVPs 


(a) function defined for all x except x = +1 


l 
l 
l 
l 
l 
l 
l 
l 
l 
l 
| fj 


ya 


(b) solution defined on interval containing x = 0 


(0, -1) 


FIGURE 1.2.4 Graphs of function 
and solution of IVP in Example 2 


INTRODUCTION TO DIFFERENTIAL EQUATIONS 


constant 3 = ce? = c. Thus y = Ze" is a solution of the IVP 


y' =y, yO) = 3. 
Now if we demand that a solution curve pass through the point (1, —2) rather than 
(0, 3), then y(1) = —2 will yield —2 = ce orc = —2e7!. In this case y= —2e* | is 
a solution of the IVP 
y =y, yl) = —2. 
The two solution curves are shown in dark blue and dark red in Figure 1.2.3. a 
The next example illustrates another first-order initial-value problem. In this 


example notice how the interval J of definition of the solution y(x) depends on the 
initial condition y(xo) = Yo. 


l EXAMPLE 2 Interval J of Definition of a Solution 


In Problem 6 of Exercises 2.2 you will be asked to show that a one-parameter family 
of solutions of the first-order differential equation vi + 2xy? = 0 is y = 1/(x? + c). 
If we impose the initial condition y(0) = —1, then substituting x = 0 and y = —1 
into the family of solutions gives —1 = 1/c or c = —1. Thus y = 1/(x? — 1). We 
now emphasize the following three distinctions: 


e Considered as a function, the domain of y = 1/(x? — 1) is the set of real 
numbers x for which y(x) is defined; this is the set of all real numbers 
except x = —1 and x = 1. See Figure 1.2.4(a). 

e Considered as a solution of the differential equation y' + 2xy? = 0, the 
interval J of definition of y = 1/(x? — 1) could be taken to be any 
interval over which y(x) is defined and differentiable. As can be seen in 
Figure 1.2.4(a), the largest intervals on which y = 1/ (x? — 1) is a solution 
are (—œ,— 1), (—1, 1), and (1, ©). 

e Considered as a solution of the initial-value problem y' + 2xy” = 0, 

y(0) = —1, the interval J of definition of y = 1/(x* — 1) could be taken to 
be any interval over which y(x) is defined, differentiable, and contains the 
initial point x = 0; the largest interval for which this is true is (—1, 1). See 
the red curve in Figure 1.2.4(b). a 


See Problems 3—6 in Exercises 1.2 for a continuation of Example 2. 


l EXAMPLE 3 Second-Order IVP 


In Example 4 of Section 1.1 we saw that x = cı cos 4t + c3 sin 4t is a two-parameter 
family of solutions of x” + 16x = 0. Find a solution of the initial-value problem 


x" + 16x = 0, A7 = -2, (3) =1. (4) 


SOLUTION We first apply x(ar/2) = —2 to the given family of solutions: cı cos 27r + 
c2 sin 277 = —2. Since cos 27 = 1 and sin 27 = 0, we find that ec = —2. We next apply 
x'(ar/2) = 1 to the one-parameter family x( = —2 cos 4t + co sin 4t. Differentiating 
and then setting t = 7/2 and x’ = 1 gives 8 sin 2a + 4c2 cos 27 = 1, from which we 
see that c) = H Hence x = —2 cos 4t + 1 sin 4ź is a solution of (4). E 


EXISTENCE AND UNIQUENESS Two fundamental questions arise in consider- 
ing an initial-value problem: 


Does a solution of the problem exist? 
Ifa solution exists, is it unique? 
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(0, 0) 


FIGURE 1.2.5 Two solutions 


of the same IVP 


(Xo, Yo) 


a kK— Io 


> 


bx 


FIGURE 1.2.6 Rectangular region R 
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For the first-order initial-value problem (2) we ask: 


Does the differential equation dy/dx = f(x, y) possess solutions? 


SE Es any of the solution curves pass through the point (xo, yo)? 


When can we be certain that there is precisely one solution curve 


Fees Bëss through the point (xo, yo)? 


Note that in Examples 1 and 3 the phrase “a solution” is used rather than “the solu- 
tion” of the problem. The indefinite article “a” is used deliberately to suggest the 
possibility that other solutions may exist. At this point it has not been demonstrated 
that there is a single solution of each problem. The next example illustrates an initial- 
value problem with two solutions. 


| EXAMPLE 4 An IVP Can Have Several Solutions 


Each of the functions y=0 and y = xx satisfies the differential equation 


dy/dx = xy!” and the initial condition y(0) = 0, so the initial-value problem 


d 
T=, yO) =0 
X 


has at least two solutions. As illustrated in Figure 1.2.5, the graphs of both functions 
pass through the same point (0, 0). D 


Within the safe confines of a formal course in differential equations one can be 
fairly confident that most differential equations will have solutions and that solutions of 
initial-value problems will probably be unique. Real life, however, is not so idyllic. 
Therefore it is desirable to know in advance of trying to solve an initial-value problem 
whether a solution exists and, when it does, whether it is the only solution of the prob- 
lem. Since we are going to consider first-order differential equations in the next two 
chapters, we state here without proof a straightforward theorem that gives conditions 
that are sufficient to guarantee the existence and uniqueness of a solution of a first-order 
initial-value problem of the form given in (2). We shall wait until Chapter 4 to address 
the question of existence and uniqueness of a second-order initial-value problem. 


THEOREM 1.2.1 Existence of a Unique Solution 


Let R be a rectangular region in the xy-plane defined bya S x S b, c S y S d 
that contains the point (xo, yo) in its interior. If f(x, y) and df/dy are continuous 
on R, then there exists some interval Jo: (xo — h, xo + h), h > 0, contained in 
[a, b], and a unique function y(x), defined on Jo, that is a solution of the initial- 
value problem (2). 


The foregoing result is one of the most popular existence and uniqueness theo- 
rems for first-order differential equations because the criteria of continuity of f(x, y) 
and df/dy are relatively easy to check. The geometry of Theorem 1.2.1 is illustrated 
in Figure 1.2.6. 


| EXAMPLE 5 Example 4 Revisited 


We saw in Example 4 that the differential equation dy/dx = xy"? possesses at least 
two solutions whose graphs pass through (0, 0). Inspection of the functions 
1/2 


F(x, y) = xy and e 
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shows that they are continuous in the upper half-plane defined by y > 0. Hence 
Theorem 1.2.1 enables us to conclude that through any point (xo, yo), yo > O in the 
upper half-plane there is some interval centered at xg on which the given differential 
equation has a unique solution. Thus, for example, even without solving it, we know 
that there exists some interval centered at 2 on which the initial-value problem 
dy/dx = xy", y(2) = 1 has a unique solution. a 


In Example 1, Theorem 1.2.1 guarantees that there are no other solutions of the 
initial-value problems vi = y, y(0) = 3 and vi = y, y(1) = —2 other than y = Zei 
and y = —2e*—!, respectively. This follows from the fact that f(x, y) = y and 
of/dy = 1 are continuous throughout the entire xy-plane. It can be further shown that 
the interval J on which each solution is defined is (—~, %). 


INTERVAL OF EXISTENCE/UNIQUENESS Suppose y(x) represents a solution 
of the initial-value problem (2). The following three sets on the real x-axis may not 
be the same: the domain of the function y(x), the interval J over which the solution 
y(x) is defined or exists, and the interval Jọ of existence and uniqueness. Example 2 
of Section 1.1 illustrated the difference between the domain of a function and the 
interval Z of definition. Now suppose (xo, yo) is a point in the interior of the rectan- 
gular region R in Theorem 1.2.1. It turns out that the continuity of the function 
f(x, y) on R by itself is sufficient to guarantee the existence of at least one solution 
of dy/dx = f(x, y), y(xo) = yo, defined on some interval 7. The interval 7 of defini- 
tion for this initial-value problem is usually taken to be the largest interval contain- 
ing x9 over which the solution y(x) is defined and differentiable. The interval / 
depends on both f(x, y) and the initial condition y(x9) = yo. See Problems 31-34 in 
Exercises 1.2. The extra condition of continuity of the first partial derivative of/dy 
on R enables us to say that not only does a solution exist on some interval Jọ con- 
taining xo, but it is the only solution satisfying y(xo) = yo. However, Theorem 1.2.1 
does not give any indication of the sizes of intervals J and Io; the interval I of 
definition need not be as wide as the region R, and the interval Ig of existence and 
uniqueness may not be as large as I. The number h > 0 that defines the interval 
Io: (xo — h, xo + h) could be very small, so it is best to think that the solution y(x) 
is unique in a local sense —that is, a solution defined near the point (xo, yo). See 
Problem 44 in Exercises 1.2. 


l REMARKS 


(i) The conditions in Theorem 1.2.1 are sufficient but not necessary. This means 
that when f(x, y) and df/dy are continuous on a rectangular region R, it must 
always follow that a solution of (2) exists and is unique whenever (xo, yo) is a 
point interior to R. However, if the conditions stated in the hypothesis of 
Theorem 1.2.1 do not hold, then anything could happen: Problem (2) may still 
have a solution and this solution may be unique, or (2) may have several solu- 
tions, or it may have no solution at all. A rereading of Example 5 reveals that the 
hypotheses of Theorem 1.2.1 do not hold on the line y = 0 for the differential 
equation dy/dx = xy", so it is not surprising, as we saw in Example 4 of this 
section, that there are two solutions defined on a common interval ~h < x < h 
satisfying y(0) = 0. On the other hand, the hypotheses of Theorem 1.2.1 do 
not hold on the line y = 1 for the differential equation dy/dx = |y — 1]. 
Nevertheless it can be proved that the solution of the initial-value problem 
dy/dx = |y — 1|, y() = 1, is unique. Can you guess this solution? 

(ii) You are encouraged to read, think about, work, and then keep in mind 
Problem 43 in Exercises 1.2. 
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EXERCISES 1.2 


Answers to selected odd-numbered problems begin on page ANS-1. 


In Problems 1 and 2, y = 1/(1 + cje *) is a one-parameter 
family of solutions of the first-order DE y' = y — wl, Find a 
solution of the first-order IVP consisting of this differential 
equation and the given initial condition. 


1. y(0) = —4 2. y(-1) =2 


In Problems 3-6, y = 1/(x” + c) is a one-parameter family 
of solutions of the first-order DE vi + 2xy* = 0. Find a 
solution of the first-order IVP consisting of this differential 
equation and the given initial condition. Give the largest 
interval J over which the solution is defined. 

3. y(2) = 3 

5. y(O) = 1 
In Problems 7-10, x = cı cos t + co sin t is a two-parameter 
family of solutions of the second-order DE x” + x = 0. Find 
a solution of the second-order IVP consisting of this differ- 
ential equation and the given initial conditions. 

7. x(0) = —1, x'(0)=8 

8. x(7/2) =0, x'(7/2)= 1 

9. x(7/6) =}, x'(m/6) =0 
10. ail = V2, x'(m/4) = 2V2 
In Problems 11-14, y = cie” + coe * is a two-parameter 
family of solutions of the second-order DE y” — y = 0. Find 
a solution of the second-order IVP consisting of this differ- 
ential equation and the given initial conditions. 
11. yO) = 1, y’(0) =2 
12. y1)=0, y’d)=e 
13. y(-1) =5, y(-l)=—-5 
14. y(0)=0, y'(0) =0 
In Problems 15 and 16 determine by inspection at least two 
solutions of the given first-order IVP. 
15. y' = 3y”, wi =0 
16. xy’ =2y, y(0)=0 
In Problems 17-24 determine a region of the xy-plane for 
which the given differential equation would have a unique 


solution whose graph passes through a point (xo, yo) in the 
region. 


19. x— = 
* dx y 


21. (4-—y?)y' = x? 
23. (x2 + y*)y' = y? 


22. (1 + y?y’ = x? 
24. (y—-x)y' =yt+x 


In Problems 25-28 determine whether Theorem 1.2.1 guar- 
antees that the differential equation y’ = Vy? — 9 pos- 
sesses a unique solution through the given point. 


25. (1, 4) 
27. (2, —3) 


26. (5,3) 
28. (—1,1) 


29. (a) By inspection find a one-parameter family of solu- 
tions of the differential equation xy’ = y. Verify that 
each member of the family is a solution of the 
initial-value problem xy’ = y, y(O) = 0. 

(b) Explain part (a) by determining a region R in the 
xy-plane for which the differential equation xy’ = y 
would have a unique solution through a point (xo, yo) 


in R. 
(c) Verify that the piecewise-defined function 
kees 
y= 
x SH 


satisfies the condition y(0) = 0. Determine whether 
this function is also a solution of the initial-value 
problem in part (a). 


30. (a) Verify that y = tan (x + c) is a one-parameter family 
of solutions of the differential equation vi = 1 + y?. 
(b) Since f(x, y) = 1 + y? and af/dy = 2y are continu- 
ous everywhere, the region R in Theorem 1.2.1 can 
be taken to be the entire xy-plane. Use the family of 
solutions in part (a) to find an explicit solution of 
the first-order initial-value problem vi = 1 + y?, 
y(0) = 0. Even though x9 = 0 is in the interval 
(—2, 2), explain why the solution is not defined on 
this interval. 


(c 


wa 


Determine the largest interval J of definition for the 
solution of the initial-value problem in part (b). 


31. (a) Verify that y= —1/(x + c) is a one-parameter 
family of solutions of the differential equation 
y a vi 
(b) Since f(x, y) = y? and df/dy = 2y are continuous 
everywhere, the region R in Theorem 1.2.1 can be 
taken to be the entire xy-plane. Find a solution from 
the family in part (a) that satisfies y(0) = 1. Then 
find a solution from the family in part (a) that 


satisfies y(0) = —1. Determine the largest interval 7 
of definition for the solution of each initial-value 
problem. 
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(c) Determine the largest interval J of definition for the 
solution of the first-order initial-value problem 
y’ = y?, y(0) = 0. [Hint: The solution is not a mem- 
ber of the family of solutions in part (a).] 


32. (a) Show that a solution from the family in part (a) 
of Problem 31 that satisfies vi = y*, y(1) = 1, is 
y=1/2-x). 

(b) Then show that a solution from the family in part (a) 
of Problem 31 that satisfies y’ = y’, y(3) = —1, is 
DEIER? 

Are the solutions in parts (a) and (b) the same? 


(c 
33. (a 


wa 


Verify that 3x? — y? = c is a one-parameter fam- 
ily of solutions of the differential equation 
y dy/dx = 3x. 

(b) By hand, sketch the graph of the implicit solution 
3x? — y? = 3. Find all explicit solutions y = $ (x) of 
the DE in part (a) defined by this relation. Give the 
interval J of definition of each explicit solution. 


Kei 


The point (—2, 3) is on the graph of 3x? — y? = 3, 
but which of the explicit solutions in part (b) satis- 
fies y(—2) = 3? 


(c 


wa 


34. (a) Use the family of solutions in part (a) of Problem 33 
to find an implicit solution of the initial-value 
problem y dy/dx = 3x, y(2) = —4. Then, by hand, 
sketch the graph of the explicit solution of this 
problem and give its interval J of definition. 


(b) Are there any explicit solutions of y dy/dx = 3x 
that pass through the origin? 


In Problems 35-38 the graph of a member of a family 
of solutions of a second-order differential equation 
d’y/dx? = f(x, y, y') is given. Match the solution curve with 
at least one pair of the following initial conditions. 

a steck y'(1)= -2 

b) y1) =0, y'=) = -—4 

© y¥)=1, vir? 

(d) yO)=—1, y'@)=2 

(e) yO) = —1, y'O)=0 

©) y0) = —4, y'(0)= -2 


35. 


5+ 


FIGURE 1.2.7 Graph for Problem 35 


36. 


37. 


38. 


FIGURE 1.2.8 Graph for Problem 36 


FIGURE 1.2.9 Graph for Problem 37 


FIGURE 1.2.10 Graph for Problem 38 


Discussion Problems 


In Problems 39 and 40 use Problem 51 in Exercises 1.1 and 
(2) and (3) of this section. 


39. 


40. 


41. 


Find a function y = f(x) whose graph at each point (x, y) 
has the slope given by 8e7*+ 6x and has the 
y-intercept (0, 9). 


Find a function y = f(x) whose second derivative is 
y” = 12x — 2 at each point (x, y) on its graph and 
y = —x + 5 is tangent to the graph at the point corre- 
sponding to x = 1. 


Consider the initial-value problem y'= x — 2y, 
y(0) = i. Determine which of the two curves shown 
in Figure 1.2.11 is the only plausible solution curve. 
Explain your reasoning. 
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43. 


44. 
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FIGURE 1.2.11 Graphs for Problem 41 


Determine a plausible value of vu for which the 
graph of the solution of the initial-value problem 
vi + 2y = 3x — 6, y(xo) = 0 is tangent to the x-axis at 
(xo, 0). Explain your reasoning. 


Suppose that the first-order differential equation 
dy/dx = f(x, y) possesses a one-parameter family of 
solutions and that f(x, y) satisfies the hypotheses of 
Theorem 1.2.1 in some rectangular region R of the 
xy-plane. Explain why two different solution curves 
cannot intersect or be tangent to each other at a point 
(xo, Yo) in R. 


The functions y(x) = jet, —% < x < % and 
0, x <0 
y(x) = Fé x=0 


have the same domain but are clearly different. See 
Figures 1.2.12(a) and 1.2.12(b), respectively. Show that 
both functions are solutions of the initial-value problem 


(a) (b) 
FIGURE 1.2.12 Two solutions of the IVP in Problem 44 


dy/dx = xy"*, y(2)=1 on the interval (—%, ©). 
Resolve the apparent contradiction between this fact 
and the last sentence in Example 5. 


Mathematical Model 


45. Population Growth Beginning in the next section 


we will see that differential equations can be used to 
describe or model many different physical systems. In 
this problem suppose that a model of the growing popu- 
lation of a small community is given by the initial-value 
problem 


dP 

a 0.15P(4) + 20, P(O) = 100, 
where P is the number of individuals in the community 
and time ¢ is measured in years. How fast—that is, at 
what rate—is the population increasing at tf = 0? How 
fast is the population increasing when the population 
is 500? 


1.3 DIFFERENTIAL EQUATIONS AS MATHEMATICAL MODELS 


REVIEW MATERIAL 


Newton’s second law of motion 
Hooke’s law 

Kirchhoff’s laws 

Archimedes’ principle 


INTRODUCTION 


these models in Chapters 3 and 5. 


MATHEMATICAL MODELS 


Units of measurement for weight, mass, and density 


In this section we introduce the notion of a differential equation as a 
mathematical model and discuss some specific models in biology, chemistry, and physics. Once we 
have studied some methods for solving DEs in Chapters 2 and 4, we return to, and solve, some of 


It is often desirable to describe the behavior of 


some real-life system or phenomenon, whether physical, sociological, or even eco- 
nomic, in mathematical terms. The mathematical description of a system of phenom- 
enon is called a mathematical model and is constructed with certain goals in mind. 
For example, we may wish to understand the mechanisms of a certain ecosystem by 
studying the growth of animal populations in that system, or we may wish to date 
fossils by analyzing the decay of a radioactive substance either in the fossil or in the 
stratum in which it was discovered. 
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Construction of a mathematical model of a system starts with 


(i) identification of the variables that are responsible for changing the 
system. We may choose not to incorporate all these variables into the 
model at first. In this step we are specifying the level of resolution of 
the model. 


Next 


(ii) we make a set of reasonable assumptions, or hypotheses, about the 
system we are trying to describe. These assumptions will also include 
any empirical laws that may be applicable to the system. 


For some purposes it may be perfectly within reason to be content with low- 
resolution models. For example, you may already be aware that in beginning 
physics courses, the retarding force of air friction is sometimes ignored in modeling 
the motion of a body falling near the surface of the Earth, but if you are a scientist 
whose job it is to accurately predict the flight path of a long-range projectile, 
you have to take into account air resistance and other factors such as the curvature 
of the Earth. 

Since the assumptions made about a system frequently involve a rate of change 
of one or more of the variables, the mathematical depiction of all these assumptions 
may be one or more equations involving derivatives. In other words, the mathemat- 
ical model may be a differential equation or a system of differential equations. 

Once we have formulated a mathematical model that is either a differential equa- 
tion or a system of differential equations, we are faced with the not insignificant 
problem of trying to solve it. Jf we can solve it, then we deem the model to be reason- 
able if its solution is consistent with either experimental data or known facts about 
the behavior of the system. But if the predictions produced by the solution are poor, 
we can either increase the level of resolution of the model or make alternative as- 
sumptions about the mechanisms for change in the system. The steps of the model- 
ing process are then repeated, as shown in the following diagram: 


- Express assumptions in terms Mathematical 
E of differential equations formulation 
l 
If necessary, 
alter assumptions Solve the DEs 
or increase resolution 
of model 
I 
l 
Check SE Display model predictions Obtain 
E (e.g., graphically) solutions 
known facts 


Of course, by increasing the resolution, we add to the complexity of the mathemati- 
cal model and increase the likelihood that we cannot obtain an explicit solution. 

A mathematical model of a physical system will often involve the variable time t. 
A solution of the model then gives the state of the system; in other words, the values 
of the dependent variable (or variables) for appropriate values of t describe the system 
in the past, present, and future. 


POPULATION DYNAMICS One of the earliest attempts to model human popula- 
tion growth by means of mathematics was by the English economist Thomas Malthus 
in 1798. Basically, the idea behind the Malthusian model is the assumption that the rate 
at which the population of a country grows at a certain time is proportional” to the total 
population of the country at that time. In other words, the more people there are at time f, 
the more there are going to be in the future. In mathematical terms, if P(t) denotes the 


“If two quantities u and v are proportional, we write u % v. This means that one quantity is a constant 
multiple of the other: u = kv. 
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total population at time ż, then this assumption can be expressed as 


dP dP 
— «P or — = kP, (1) 
dt dt 

where k is a constant of proportionality. This simple model, which fails to take into 
account many factors that can influence human populations to either grow or decline 
(immigration and emigration, for example), nevertheless turned out to be fairly accu- 
rate in predicting the population of the United States during the years 1790—1860. 
Populations that grow at a rate described by (1) are rare; nevertheless, (1) is still used 
to model growth of small populations over short intervals of time (bacteria growing 
in a petri dish, for example). 


RADIOACTIVE DECAY The nucleus of an atom consists of combinations of pro- 
tons and neutrons. Many of these combinations of protons and neutrons are unstable — 
that is, the atoms decay or transmute into atoms of another substance. Such nuclei are 
said to be radioactive. For example, over time the highly radioactive radium, Ra-226, 
transmutes into the radioactive gas radon, Rn-222. To model the phenomenon of 
radioactive decay, it is assumed that the rate dA/dt at which the nuclei of a sub- 
stance decay is proportional to the amount (more precisely, the number of nuclei) 
A(t) of the substance remaining at time t: 


dA _ 


dA 
a A — 
dt 


J kA. (2) 
Of course, equations (1) and (2) are exactly the same; the difference is only in the in- 
terpretation of the symbols and the constants of proportionality. For growth, as we 
expect in (1), k > 0, and for decay, as in (2), k < 0. 

The model (1) for growth can also be seen as the equation dS/dt = rS, which 
describes the growth of capital 5 when an annual rate of interest r is compounded 
continuously. The model (2) for decay also occurs in biological applications such as 
determining the half-life of a drug—the time that it takes for 50% of a drug to be 
eliminated from a body by excretion or metabolism. In chemistry the decay model 
(2) appears in the mathematical description of a first-order chemical reaction. The 
point is this: 


A single differential equation can serve as a mathematical model for many 
different phenomena. 


Mathematical models are often accompanied by certain side conditions. For ex- 
ample, in (1) and (2) we would expect to know, in turn, the initial population Po and 
the initial amount of radioactive substance Ag on hand. If the initial point in time is 
taken to be t = 0, then we know that P(0) = Po and A(0) = Ao. In other words, a 
mathematical model can consist of either an initial-value problem or, as we shall see 
later on in Section 5.2, a boundary-value problem. 


NEWTON’S LAW OF COOLING/WARMING According to Newton’s empiri- 
cal law of cooling/warming, the rate at which the temperature of a body changes is 
proportional to the difference between the temperature of the body and the temper- 
ature of the surrounding medium, the so-called ambient temperature. If T(t) repre- 
sents the temperature of a body at time t, T„ the temperature of the surrounding 
medium, and dT/dt the rate at which the temperature of the body changes, then 
Newton’s law of cooling/warming translates into the mathematical statement 


dT dT KT — T,) (3) 
SN or Eege = $ 
dt m dt m 


where k is a constant of proportionality. In either case, cooling or warming, if Tm is a 
constant, it stands to reason that k < 0. 
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SPREAD OF A DISEASE A contagious disease—for example, a flu virus—is 
spread throughout a community by people coming into contact with other people. Let 
x(t) denote the number of people who have contracted the disease and y(f) denote the 
number of people who have not yet been exposed. It seems reasonable to assume that 
the rate dx/dt at which the disease spreads is proportional to the number of encoun- 
ters, or interactions, between these two groups of people. If we assume that the num- 
ber of interactions is jointly proportional to x(t) and y(t)— that is, proportional to the 
product xy—then 


— = Kr, (4) 


where k is the usual constant of proportionality. Suppose a small community has a 
fixed population of n people. If one infected person is introduced into this commu- 
nity, then it could be argued that x(t) and vi are related by x + y =n + 1. Using 
this last equation to eliminate y in (4) gives us the model 


d 
= = kein + 1 = xy. (5) 


An obvious initial condition accompanying equation (5) is x(0) = 1. 


CHEMICAL REACTIONS The disintegration of a radioactive substance, governed 
by the differential equation (1), is said to be a first-order reaction. In chemistry 
a few reactions follow this same empirical law: If the molecules of substance A 
decompose into smaller molecules, it is a natural assumption that the rate at which 
this decomposition takes place is proportional to the amount of the first substance 
that has not undergone conversion; that is, if X(f) is the amount of substance A 
remaining at any time, then dX/dt = kX, where k is a negative constant since X is 
decreasing. An example of a first-order chemical reaction is the conversion of t-butyl 
chloride, (CH: SCC, into t-butyl alcohol, (CH3)3COH: 


(CH;);CCl + NaOH — (CH;);COH + NaCl. 


Only the concentration of the t-butyl chloride controls the rate of reaction. But in the 
reaction 


CHOCO + NaOH — CHOH + NaCl 


one molecule of sodium hydroxide, NaOH, is consumed for every molecule of 
methyl chloride, CH3Cl, thus forming one molecule of methyl alcohol, CH3OH, and 
one molecule of sodium chloride, NaCl. In this case the rate at which the reaction 
proceeds is proportional to the product of the remaining concentrations of CH3Cl and 
NaOH. To describe this second reaction in general, let us suppose one molecule of a 
substance A combines with one molecule of a substance B to form one molecule of a 
substance C. If X denotes the amount of chemical C formed at time ¢ and if œ and B 
are, in turn, the amounts of the two chemicals A and B at t = O (the initial amounts), 
then the instantaneous amounts of A and B not converted to chemical C are a — X 
and B — X, respectively. Hence the rate of formation of C is given by 


dX 

Pra k(a — XXP — X), (6) 
t 

where k is a constant of proportionality. A reaction whose model is equation (6) is 

said to be a second-order reaction. 


MIXTURES The mixing of two salt solutions of differing concentrations gives 
rise to a first-order differential equation for the amount of salt contained in the mix- 
ture. Let us suppose that a large mixing tank initially holds 300 gallons of brine (that 
is, water in which a certain number of pounds of salt has been dissolved). Another 
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input rate of brine 
3 gal/min 


constant 
300 gal 


output rate of brine 
3 gal/min 


FIGURE 1.3.1 Mixing tank 


FIGURE 1.3.2 Draining tank 
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brine solution is pumped into the large tank at a rate of 3 gallons per minute; the 
concentration of the salt in this inflow is 2 pounds per gallon. When the solution in 
the tank is well stirred, it is pumped out at the same rate as the entering solution. See 
Figure 1.3.1. If A(t) denotes the amount of salt (measured in pounds) in the tank at 
time ż, then the rate at which A(t) changes is a net rate: 


dA input rate output rate 
= ~ = Rin ge Rou (7) 
dt of salt of salt 


The input rate Rj, at which salt enters the tank is the product of the inflow concentra- 
tion of salt and the inflow rate of fluid. Note that Rin is measured in pounds per 
minute: 


concentration 
of salt input rate input rate 
in inflow of brine of salt 


R,, = (2 Ib/gal) - (3 gal/min) = (6 Ib/min). 


Now, since the solution is being pumped out of the tank at the same rate that it is 
pumped in, the number of gallons of brine in the tank at time ¢ is a constant 300 gal- 
lons. Hence the concentration of the salt in the tank as well as in the outflow is 
c(t) = A(t)/300 Ib/gal, so the output rate R ou of salt is 


concentration 
of salt output rate output rate 
in outflow of brine of salt 
l | l 
_ (AM) : .._ A(t) ; 
Row = E sel (3 gal/min) = T00 Ib/min. 
The net rate (7) then becomes 
dA A dA | 
Én ke hn A -F =; A = 6. 8 
di 100 Ta 100 Ge 


If Fin and Four denote general input and output rates of the brine solutions,’ then 
there are three possibilities: rin = Fouts Fin > Your, ANd Fin < ou In the analysis lead- 
ing to (8) we have assumed that rj, = Fout. In the latter two cases the number of gal- 
lons of brine in the tank is either increasing (Fin > Fout) or decreasing (Fin < Four) at 
the net rate rin — Fout. See Problems 10—12 in Exercises 1.3. 


DRAINING A TANK In hydrodynamics Torricelli’s law states that the speed v of 
efflux of water though a sharp-edged hole at the bottom of a tank filled to a depth h 
is the same as the speed that a body (in this case a drop of water) would acquire in 
falling freely from a height h—that is, v = V2gh, where g is the acceleration due to 
gravity. This last expression comes from equating the kinetic energy sm v? with the 
potential energy mgh and solving for v. Suppose a tank filled with water is allowed to 
drain through a hole under the influence of gravity. We would like to find the depth h 
of water remaining in the tank at time t. Consider the tank shown in Figure 1.3.2. If 
the area of the hole is A, (in ft?) and the speed of the water leaving the tank is 
v = V2gh (in ft/s), then the volume of water leaving the tank per second is A, V2gh 
(in ft?/s). Thus if V(t) denotes the volume of water in the tank at time f, then 


dV 


ae EE (9) 


“Don’t confuse these symbols with Ra, and Rou, which are input and output rates of salt. 
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o 
C 


(a) LRC-series circuit 


Inductor 

inductance L: henries (h) 
di 
dt 


L 


voltage drop across: L 


Resistor 
resistance R: ohms (Q) 
voltage drop across: iR 


Capacitor 
capacitance C: farads (f) 


voltage drop across: cd 


| 

| 

C 
(b) 

FIGURE 1.3.3 Symbols, units, and 

voltages. Current i(t) and charge q(t) are 


measured in amperes (A) and coulombs 
(C), respectively 


FIGURE 1.3.4 Position of rock 


measured from ground level 


INTRODUCTION TO DIFFERENTIAL EQUATIONS 


where the minus sign indicates that V is decreasing. Note here that we are ignoring 
the possibility of friction at the hole that might cause a reduction of the rate of flow 
there. Now if the tank is such that the volume of water in it at time ¢ can be written 
V(t) = Ayh, where A, (in ft?) is the constant area of the upper surface of the water 
(see Figure 1.3.2), then dV/dt = A, dh/dt. Substituting this last expression into (9) 
gives us the desired differential equation for the height of the water at time f: 


2S 


2gh. 10 
gë A, 8 (10) 


It is interesting to note that (10) remains valid even when A, is not constant. In this 
case we must express the upper surface area of the water as a function of h— that is, 
Aw = A(h). See Problem 14 in Exercises 1.3. 


SERIES CIRCUITS Consider the single-loop series circuit shown in Figure 1.3.3(a), 
containing an inductor, resistor, and capacitor. The current in a circuit after a switch 
is closed is denoted by i(t); the charge on a capacitor at time t is denoted by q(t). The 
letters L, R, and C are known as inductance, resistance, and capacitance, respectively, 
and are generally constants. Now according to Kirchhoff’s second law, the im- 
pressed voltage E(f) on a closed loop must equal the sum of the voltage drops in the 
loop. Figure 1.3.3(b) shows the symbols and the formulas for the respective voltage 
drops across an inductor, a capacitor, and a resistor. Since current i(f) is related to 
charge q(t) on the capacitor by i = dq/dt, adding the three voltages 


inductor resistor capacitor 
d d? d 1 
LŽ” =L Baba 44 
dt dr dr C 


and equating the sum to the impressed voltage yields a second-order differential 
equation 


d’q doa. 
L ar gia Es = E(t). 11 
TE aa 4 (t) (11) 


We will examine a differential equation analogous to (11) in great detail in 
Section 5.1. 


FALLING BODIES To construct a mathematical model of the motion of a body 
moving in a force field, one often starts with Newton’s second law of motion. Recall 
from elementary physics that Newton’s first law of motion states that a body either 
will remain at rest or will continue to move with a constant velocity unless acted on 
by an external force. In each case this is equivalent to saying that when the sum of 
the forces F = X, F,—that is, the net or resultant force—acting on the body is zero, 
then the acceleration a of the body is zero. Newton’s second law of motion 
indicates that when the net force acting on a body is not zero, then the net force is 
proportional to its acceleration a or, more precisely, F = ma, where m is the mass of 
the body. 

Now suppose a rock is tossed upward from the roof of a building as illustrated in 
Figure 1.3.4. What is the position st) of the rock relative to the ground at time t? The 
acceleration of the rock is the second derivative d*s /dt”. If we assume that the up- 
ward direction is positive and that no force acts on the rock other than the force of 
gravity, then Newton’s second law gives 


d? d' 
jt = —mg or T = —g. (12) 
In other words, the net force is simply the weight F = Fı = —W of the rock near the 
surface of the Earth. Recall that the magnitude of the weight is W = mg, where m is 
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kv 
positive air resistance 
direction @ 
| gravity 


mg 
FIGURE 1.3.5 Falling body of mass m 


(a) suspension bridge cable 


(b) telephone wires 


FIGURE 1.3.6 Cables suspended 
between vertical supports 


FIGURE 1.3.7 Element of cable 
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the mass of the body and g is the acceleration due to gravity. The minus sign in (12) is 
used because the weight of the rock is a force directed downward, which is opposite 
to the positive direction. If the height of the building is sọ and the initial velocity of the 
rock is vo, then s is determined from the second-order initial-value problem 


es g SO) = so -s (0) = vo: (13) 


Although we have not been stressing solutions of the equations we have con- 
structed, note that (13) can be solved by integrating the constant —g twice with 
respect to t The initial conditions determine the two constants of integration. 
From elementary physics you might recognize the solution of (13) as the formula 
s(t) = —3gt? + vot + so. 


FALLING BODIES AND AIR RESISTANCE Before Galileo’s famous experiment 
from the leaning tower of Pisa, it was generally believed that heavier objects in free 
fall, such as a cannonball, fell with a greater acceleration than lighter objects, such as 
a feather. Obviously, a cannonball and a feather when dropped simultaneously from 
the same height do fall at different rates, but it is not because a cannonball is heavier. 
The difference in rates is due to air resistance. The resistive force of air was ignored 
in the model given in (13). Under some circumstances a falling body of mass m, such 
as a feather with low density and irregular shape, encounters air resistance propor- 
tional to its instantaneous velocity v. If we take, in this circumstance, the positive 
direction to be oriented downward, then the net force acting on the mass is given by 
F = F; + F, = mg — kv, where the weight F = mg of the body is force acting in the 
positive direction and air resistance F) = —kv is a force, called viscous damping, 
acting in the opposite or upward direction. See Figure 1.3.5. Now since v is related to 
acceleration a by a = dv/dt, Newton’s second law becomes F = ma = m dv/dt. By 
equating the net force to this form of Newton’s second law, we obtain a first-order 
differential equation for the velocity v(t) of the body at time t, 


Ep (14) 
we ns i 


Here k is a positive constant of proportionality. If s(t) is the distance the body falls in 
time £ from its initial point of release, then v = ds /dt and a = dv/dt = d'eidri. In 
terms of s, (14) is a second-order differential equation 


maa "s -ky or nmay ea e (15) 


SUSPENDED CABLES Suppose a flexible cable, wire, or heavy rope is suspended 
between two vertical supports. Physical examples of this could be one of the two 
cables supporting the roadbed of a suspension bridge as shown in Figure 1.3.6(a) or 
a long telephone wire strung between two posts as shown in Figure 1.3.6(b). Our goal 
is to construct a mathematical model that describes the shape that such a cable 
assumes. 

To begin, let’s agree to examine only a portion or element of the cable between 
its lowest point P; and any arbitrary point P2. As drawn in blue in Figure 1.3.7, this 
element of the cable is the curve in a rectangular coordinate system with y-axis cho- 
sen to pass through the lowest point P; on the curve and the x-axis chosen a units 
below P. Three forces are acting on the cable: the tensions T; and T; in the cable 
that are tangent to the cable at P; and P2, respectively, and the portion W of the total 
vertical load between the points P and P Let Ti =|T,|, Tə = |T|, and 
W = |W] denote the magnitudes of these vectors. Now the tension TT: resolves 
into horizontal and vertical components (scalar quantities) T cos 0 and T sin 0. 
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Because of static equilibrium we can write 
T, = T, cos 0 and W = T, sin 0. 


By dividing the last equation by the first, we eliminate T, and get tan 8 = W/T. But 
because dy/dx = tan 0, we arrive at 


dy W 
—= = 16 
de Ti (16) 
This simple first-order differential equation serves as a model for both the shape of a 
flexible wire such as a telephone wire hanging under its own weight and the shape of 
the cables that support the roadbed of a suspension bridge. We will come back to 
equation (16) in Exercises 2.2 and Section 5.3. 


WHAT LIES AHEAD Throughout this text you will see three different types of 
approaches to, or analyses of, differential equations. Over the centuries differential 
equations would often spring from the efforts of a scientist or engineer to describe 
some physical phenomenon or to translate an empirical or experimental law into 
mathematical terms. As a consequence a scientist, engineer, or mathematician would 
often spend many years of his or her life trying to find the solutions of a DE. With a 
solution in hand, the study of its properties then followed. This quest for solutions is 
called by some the analytical approach to differential equations. Once they realized 
that explicit solutions are at best difficult to obtain and at worst impossible to obtain, 
mathematicians learned that a differential equation itself could be a font of valuable 
information. It is possible, in some instances, to glean directly from the differential 
equation answers to questions such as Does the DE actually have solutions? If a 
solution of the DE exists and satisfies an initial condition, is it the only such solu- 
tion? What are some of the properties of the unknown solutions? What can we say 
about the geometry of the solution curves? Such an approach is qualitative analysis. 
Finally, if a differential equation cannot be solved by analytical methods, yet we 
can prove that a solution exists, the next logical query is Can we somehow approxi- 
mate the values of an unknown solution? Here we enter the realm of numerical 
analysis. An affirmative answer to the last question stems from the fact that a differ- 
ential equation can be used as a cornerstone for constructing very accurate approxi- 
mation algorithms. In Chapter 2 we start with qualitative considerations of first-order 
ODEs, then examine analytical stratagems for solving some special first-order equa- 
tions, and conclude with an introduction to an elementary numerical method. See 
Figure 1.3.8. 


(a) analytical 


(b) qualitative (c) numerical 


FIGURE 1.3.8 Different approaches to the study of differential equations 
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REMARKS 


Each example in this section has described a dynamical system —a system that 
changes or evolves with the flow of time ¢. Since the study of dynamical 
systems is a branch of mathematics currently in vogue, we shall occasionally 
relate the terminology of that field to the discussion at hand. 

In more precise terms, a dynamical system consists of a set of time- 
dependent variables, called state variables, together with a rule that enables 
us to determine (without ambiguity) the state of the system (this may be a past, 
present, or future state) in terms of a state prescribed at some time tọ. Dynamical 
systems are classified as either discrete-time systems or continuous-time systems. 
In this course we shall be concerned only with continuous-time systems — 
systems in which all variables are defined over a continuous range of time. The 
rule, or mathematical model, in a continuous-time dynamical system is a differ- 
ential equation or a system of differential equations. The state of the system 
at a time f is the value of the state variables at that time; the specified state of 
the system at a time tọ is simply the initial conditions that accompany the math- 
ematical model. The solution of the initial-value problem is referred to as the 
response of the system. For example, in the case of radioactive decay, the rule 
is dA/dt = kA. Now if the quantity of a radioactive substance at some time fo is 
known, say A(tfo) = Ao, then by solving the rule we find that the response of the 
system for t = fo is A(t) = Age“) (see Section 3.1). The response A(f) is the 
single state variable for this system. In the case of the rock tossed from the roof 
of a building, the response of the system—the solution of the differential 
equation d7s/dt* = —g, subject to the initial state s(0) = so, s'(0) = vo, is the 
function s(t) = —igr + Vot + So,0 =t=T, where T represents the time 
when the rock hits the ground. The state variables are s(t) and s’(t), which 
are the vertical position of the rock above ground and its velocity at time t, 
respectively. The acceleration s”(f) is not a state variable, since we have to know 
only any initial position and initial velocity at a time fo to uniquely determine 
the rock’s position s(f) and velocity si) = vU) for any time in the interval 
to S t ST. The acceleration su) = a(t) is, of course, given by the differential 
equation su) = —g,0<t<T. 

One last point: Not every system studied in this text is a dynamical system. 
We shall also examine some static systems in which the model is a differential 
equation. 


| | EXERCISES 1.3 


Answers to selected odd-numbered problems begin on page ANS-1. 


Population Dynamics 


1. Under the same assumptions that underlie the model in 
(1), determine a differential equation for the population 
P(t) of a country when individuals are allowed to 
immigrate into the country at a constant rate r > 0. 
What is the differential equation for the population P(t) 
of the country when individuals are allowed to emigrate 
from the country at a constant rate r > 0? 


. The population model given in (1) fails to take death 
into consideration; the growth rate equals the birth rate. 
In another model of a changing population of a commu- 
nity it is assumed that the rate at which the population 
changes is a net rate—that is, the difference between 


the rate of births and the rate of deaths in the commu- 
nity. Determine a model for the population P(t) if both 
the birth rate and the death rate are proportional to the 
population present at time t. 


. Using the concept of net rate introduced in Problem 2, 


determine a model for a population P if the birth rate 
is proportional to the population present at time rf but the 
death rate is proportional to the square of the population 
present at time t. 


. Modify the model in Problem 3 for net rate at which 


the population P(t) of a certain kind of fish changes by 
also assuming that the fish are harvested at a constant 
rateh > 0. 
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Newton’s Law of Cooling/Warming 


5. 


A cup of coffee cools according to Newton’s law of 
cooling (3). Use data from the graph of the temperature 
T(t) in Figure 1.3.9 to estimate the constants Tm, To, and 
k in a model of the form of a first-order initial-value 
problem: dT/dt = k(T — Tm), T(O) = To. 


EE 


0 25 50 75 100 
minutes 


FIGURE 1.3.9 Cooling curve in Problem 5 


. The ambient temperature 7,, in (3) could be a function 


of time ¢ Suppose that in an artificially controlled 
environment, Tat?) is periodic with a 24-hour period, 
as illustrated in Figure 1.3.10. Devise a mathematical 
model for the temperature T(t) of a body within this 
environment. 


ACS F 
120 + 


1004 
804 
605 


t > 


0 12 24 36 48 t 
midnight noon midnight noon midnight 


FIGURE 1.3.10 Ambient temperature in Problem 6 


Spread of a Disease/Technology 


7. 


Suppose a student carrying a flu virus returns to an iso- 
lated college campus of 1000 students. Determine a dif- 
ferential equation for the number of people x(t) who have 
contracted the flu if the rate at which the disease spreads 
is proportional to the number of interactions between the 
number of students who have the flu and the number of 
students who have not yet been exposed to it. 


At a time denoted as ¢ = 0 a technological innovation is 
introduced into a community that has a fixed population 
of n people. Determine a differential equation for the 


INTRODUCTION TO DIFFERENTIAL EQUATIONS 


number of people x(t) who have adopted the innovation 
at time ¢ if it is assumed that the rate at which the innova- 
tions spread through the community is jointly propor- 
tional to the number of people who have adopted it and 
the number of people who have not adopted it. 


Mixtures 


9. 


10. 


11. 


12. 


Suppose that a large mixing tank initially holds 300 gal- 
lons of water in which 50 pounds of salt have been dis- 
solved. Pure water is pumped into the tank at a rate of 
3 gal/min, and when the solution is well stirred, it is 
then pumped out at the same rate. Determine a differen- 
tial equation for the amount of salt A(f) in the tank at 
time t. What is A(O)? 


Suppose that a large mixing tank initially holds 300 gal- 
lons of water is which 50 pounds of salt have been 
dissolved. Another brine solution is pumped into the tank 
at a rate of 3 gal/min, and when the solution is well 
stirred, it is then pumped out at a slower rate of 2 gal/min. 
If the concentration of the solution entering is 2 lb/gal, 
determine a differential equation for the amount of salt 
A(t) in the tank at time t. 


What is the differential equation in Problem 10, if the 
well-stirred solution is pumped out at a faster rate of 
3.5 gal/min? 


Generalize the model given in equation (8) on page 23 
by assuming that the large tank initially contains No 
number of gallons of brine, rj, and Four are the input and 
output rates of the brine, respectively (measured in gal- 
lons per minute), cj, is the concentration of the salt in 
the inflow, c(t) the concentration of the salt in the tank 
as well as in the outflow at time t (measured in pounds 
of salt per gallon), and A(t) is the amount of salt in the 
tank at time t. 


Draining a Tank 


13. 


Suppose water is leaking from a tank through a circular 
hole of area Au at its bottom. When water leaks through a 
hole, friction and contraction of the stream near the hole 
reduce the volume of water leaving the tank per second to 
cA, V2gh, where c (0 < c < 1) is an empirical constant. 
Determine a differential equation for the height / of water 
at time £ for the cubical tank shown in Figure 1.3.11. The 
radius of the hole is 2 in., and g = 32 ft/s”. 


circular =| 
hole 


FIGURE 1.3.11 Cubical tank in Problem 13 
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14. The right-circular conical tank shown in Figure 1.3.12 
loses water out of a circular hole at its bottom. Determine 
a differential equation for the height of the water h at 
time ¢. The radius of the hole is 2 in., g = 32 ft/s?, and 
the friction/contraction factor introduced in Problem 13 
isc = 0.6. 


FIGURE 1.3.12 Conical tank in Problem 14 


Series Circuits 


15. A series circuit contains a resistor and an inductor as 
shown in Figure 1.3.13. Determine a differential equa- 
tion for the current i(f) if the resistance is R, the induc- 
tance is L, and the impressed voltage is E(t). 


© 


FIGURE 1.3.13 LR series circuit in Problem 15 


16. A series circuit contains a resistor and a capacitor as 
shown in Figure 1.3.14. Determine a differential equa- 
tion for the charge q(t) on the capacitor if the resis- 
tance is R, the capacitance is C, and the impressed 
voltage is E(t). 


FIGURE 1.3.14 RC series circuit in Problem 16 


Falling Bodies and Air Resistance 


17. For high-speed motion through the air—such as the 
skydiver shown in Figure 1.3.15, falling before the para- 
chute is opened —air resistance is closer to a power of 
the instantaneous velocity v(t). Determine a differential 
equation for the velocity v(t) of a falling body of mass m 
if air resistance is proportional to the square of the 
instantaneous velocity. 


mg 


FIGURE 1.3.15 Air resistance proportional to square of 
velocity in Problem 17 


Newton’s Second Law and Archimedes’ Principle 


18. A cylindrical barrel s feet in diameter of weight w Ib 
is floating in water as shown in Figure 1.3.16(a). After 
an initial depression the barrel exhibits an up-and- 
down bobbing motion along a vertical line. Using 
Figure 1.3.16(b), determine a differential equation for 
the vertical displacement y(¢) if the origin is taken to be 
on the vertical axis at the surface of the water when the 
barrel is at rest. Use Archimedes’ principle: Buoyancy, 
or upward force of the water on the barrel, is equal to 
the weight of the water displaced. Assume that the 
downward direction is positive, that the weight density 
of water is 62.4 lb/ft?, and that there is no resistance 
between the barrel and the water. 


surface vo > YO 


(a) (b) 


FIGURE 1.3.16 Bobbing motion of floating barrel in 
Problem 18 


Newton’s Second Law and Hooke’s Law 


19. After a mass m is attached to a spring, it stretches it 
s units and then hangs at rest in the equilibrium position 
as shown in Figure 1.3.17(b). After the spring/mass 


unstretched 8 5 Geh 
spring +- ee 
equilibrium x() >0 
position | 
(a) (b) (©) 


FIGURE 1.3.17  Spring/mass system in Problem 19 
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system has been set in motion, let x(t) denote the di- 
rected distance of the mass beyond the equilibrium po- 
sition. As indicated in Figure 1.3.17(c), assume that the 
downward direction is positive, that the motion takes 
place in a vertical straight line through the center of 
gravity of the mass, and that the only forces acting on 
the system are the weight of the mass and the restoring 
force of the stretched spring. Use Hooke’s law: The 
restoring force of a spring is proportional to its total 
elongation. Determine a differential equation for the 
displacement x(t) at time t. 


In Problem 19, what is a differential equation for the 
displacement x(t) if the motion takes place in a medium 
that imparts a damping force on the spring/mass system 
that is proportional to the instantaneous velocity of the 
mass and acts in a direction opposite to that of motion? 


Newton’s Second Law and the Law 
of Universal Gravitation 


21. 


22. 


By Newton’s universal law of gravitation the free-fall 
acceleration a of a body, such as the satellite shown in 
Figure 1.3.18, falling a great distance to the surface is not 
the constant g. Rather, the acceleration a is inversely pro- 
portional to the square of the distance from the center of 
the Earth, a = k/ r?, where k is the constant of proportion- 
ality. Use the fact that at the surface of the Earth r = R and 
a = g to determine k. If the positive direction is upward, 
use Newton’s second law and his universal law of gravita- 
tion to find a differential equation for the distance z: 


satellite of 
mass m Be, 


FIGURE 1.3.18 Satellite 
in Problem 21 


Earth of mass M 


Suppose a hole is drilled through the center of the Earth 
and a bowling ball of mass m is dropped into the hole, as 
shown in Figure 1.3.19. Construct a mathematical model 
that describes the motion of the ball. At time r let r de- 
note the distance from the center of the Earth to the mass 
m, M denote the mass of the Earth, M, denote the mass of 
that portion of the Earth within a sphere of radius r, and 
6 denote the constant density of the Earth. 


surface 


FIGURE 1.3.19 Hole through 
Earth in Problem 22 
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Additional Mathematical Models 


23. 


24. 


25. 


26. 


27. 


Learning Theory In the theory of learning, the rate at 
which a subject is memorized is assumed to be pro- 
portional to the amount that is left to be memorized. 
Suppose M denotes the total amount of a subject to be 
memorized and A(f) is the amount memorized in time t. 
Determine a differential equation for the amount A(t). 


Forgetfulness In Problem 23 assume that the rate at 
which material is forgotten is proportional to the amount 
memorized in time t. Determine a differential equation 
for the amount A(t) when forgetfulness is taken into 
account. 


Infusion of a Drug A drug is infused into a patient’s 
bloodstream at a constant rate of r grams per second. 
Simultaneously, the drug is removed at a rate proportional 
to the amount x(t) of the drug present at time t. Determine 
a differential equation for the amount x(t). 


Tractrix A person P, starting at the origin, moves in the 
direction of the positive x-axis, pulling a weight along 
the curve C, called a tractrix, as shown in Figure 1.3.20. 
The weight, initially located on the y-axis at (0, s), is 
pulled by a rope of constant length s, which is kept taut 
throughout the motion. Determine a differential equation 
for the path C of motion. Assume that the rope is always 
tangent to C. 


FIGURE 1.3.20 ` Tractrix curve in Problem 26 


Reflecting Surface Assume that when the plane 
curve C shown in Figure 1.3.21 is revolved about the 
X-axis, it generates a surface of revolution with the prop- 
erty that all light rays L parallel to the x-axis striking the 
surface are reflected to a single point O (the origin). Use 
the fact that the angle of incidence is equal to the angle 
of reflection to determine a differential equation that 


tangent 


A 


> 
x 


FIGURE 1.3.21 Reflecting surface in Problem 27 


Copyright 2009 Cengage Learning, Inc. All Rights Reserved. May not be copied, scanned, or duplicated, in whole or in part. 


Licensed to: iChapters User 


1.3 DIFFERENTIAL EQUATIONS AS MATHEMATICAL MODELS ° 31 


describes the shape of the curve C. Such a curve C is 
important in applications ranging from construction of 
telescopes to satellite antennas, automobile headlights, 
and solar collectors. [Hint: Inspection of the figure 
shows that we can write dh = 20. Why? Now use an 
appropriate trigonometric identity. ] 


Discussion Problems 


28. 


29. 


30. 


31. 


32. 


Reread Problem 41 in Exercises 1.1 and then give an 
explicit solution P(t) for equation (1). Find a one- 
parameter family of solutions of (1). 


Reread the sentence following equation (3) and assume 
that T, is a positive constant. Discuss why we would ex- 
pect k < 0 in (3) in both cases of cooling and warming. 
You might start by interpreting, say, T(t) > Tm in a 
graphical manner. 


Reread the discussion leading up to equation (8). If we 
assume that initially the tank holds, say, 50 Ib of salt, it 
stands to reason that because salt is being added to the 
tank continuously for t > 0, A(t) should be an increas- 
ing function. Discuss how you might determine from 
the DE, without actually solving it, the number of 
pounds of salt in the tank after a long period of time. 


Population Model The differential equation 
dP 
E = (kcos t)P, where k is a positive constant, is a 


model of human population P(t) of a certain commu- 
nity. Discuss an interpretation for the solution of this 
equation. In other words, what kind of population do 
you think the differential equation describes? 


Rotating Fluid As shown in Figure 1.3.22(a), a right- 
circular cylinder partially filled with fluid is rotated 
with a constant angular velocity w about a vertical y-axis 
through its center. The rotating fluid forms a surface of 
revolution S. To identify S, we first establish a coordinate 
system consisting of a vertical plane determined by the 
y-axis and an x-axis drawn perpendicular to the y-axis 
such that the point of intersection of the axes (the origin) 
is located at the lowest point on the surface S. We then 
seek a function y = f(x) that represents the curve C of in- 
tersection of the surface A and the vertical coordinate 
plane. Let the point P(x, y) denote the position of a parti- 
cle of the rotating fluid of mass m in the coordinate 
plane. See Figure 1.3.22(b). 


(a) At P there is a reaction force of magnitude F due to 
the other particles of the fluid which is normal to the 
surface A. By Newton’s second law the magnitude 
of the net force acting on the particle is mw*x. What 
is this force? Use Figure 1.3.22(b) to discuss the na- 


ture and origin of the equations 
F cos 0 = mg, F sin 0 = max. 


(b) Use part (a) to find a first-order differential equation 
that defines the function y = f(x). 


(a) 


curve C of 

intersection 
of xy-plane 
and surface A 
of revolution 


ay 


tangent line to 
curve C at P 


(b) 
FIGURE 1.3.22 Rotating fluid in Problem 32 


33. Falling Body In Problem 21, suppose r= R + s, 


34. 


where s is the distance from the surface of the Earth to 
the falling body. What does the differential equation 
obtained in Problem 21 become when s is very small in 
comparison to R? [Hint: Think binomial series for 


(R+s)?7=R*(14+5/R) il 


Raindrops Keep Falling In meteorology the term 
virga refers to falling raindrops or ice particles that 
evaporate before they reach the ground. Assume that a 
typical raindrop is spherical. Starting at some time, 
which we can designate as t = 0, the raindrop of radius 
ro falls from rest from a cloud and begins to evaporate. 


(a) If it is assumed that a raindrop evaporates in such a 
manner that its shape remains spherical, then it also 
makes sense to assume that the rate at which the rain- 
drop evaporates—that is, the rate at which it loses 
mass—is proportional to its surface area. Show that 
this latter assumption implies that the rate at which 
the radius r of the raindrop decreases is a constant. 
Find r(t). [Hint: See Problem 51 in Exercises 1.1.] 


(b) If the positive direction is downward, construct a 
mathematical model for the velocity v of the falling 
raindrop at time ¢. Ignore air resistance. [Hint: 
When the mass m of an object is changing with 


d 
time, Newton’s second law becomes F = E (my), 


where F is the net force acting on the body and mv 
is its momentum. | 
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35. Let It Snow The “snowplow problem” is a classic and 
appears in many differential equations texts but was 
probably made famous by Ralph Palmer Agnew: 


“One day it started snowing at a heavy and steady 
rate. A snowplow started out at noon, going 2 miles 
the first hour and 1 mile the second hour. What time 
did it start snowing?” 


INTRODUCTION TO DIFFERENTIAL EQUATIONS 


Find the text Differential Equations, Ralph Palmer Agnew, 
McGraw-Hill Book Co., and then discuss the construction 
and solution of the mathematical model. 


36. Reread this section and classify each mathematical 
model as linear or nonlinear. 


CHAPTER 1 IN REVIEW 


Answers to selected odd-numbered problems begin on page ANS-1. 


In Problems 1 and 2 fill in the blank and then write this result 
as a linear first-order differential equation that is free of the 
symbol cı and has the form dy/dx = f(x, y). The symbol cı 
represents a constant. 


d eg 
2: ae + ce“) = 


In Problems 3 and 4 fill in the blank and then write this result 
as a linear second-order differential equation that is free of 
the symbols cı and cz and has the form F(y, y”) = 0. The 
symbols c1, c2, and k represent constants. 


2 
dx? 
2 


d 
4. ES (cı cosh kx + c, sinh kx) = 


3 


(cı cos kx + c sin kx) = 


In Problems 5 and 6 compute y’ and y” and then combine 
these derivatives with y as a linear second-order differential 
equation that is free of the symbols cı and ec: and has the form 


F(y, vi y") = 0. The symbols cı and c2 represent constants. 
5. y = cye* + coxe” 6. y = ce” cos x + coe* sin x 


In Problems 7—12 match each of the given differential equa- 
tions with one or more of these solutions: 


@)y=0, @)y=2, ©y=2x, @ y= ie 
7. xy’ = 2y 8. y' =2 

9, y =2y—-4 10. xy’ = y 

11. vi + 9y = 18 12. xy"—y' =0 


In Problems 13 and 14 determine by inspection at least one 
solution of the given differential equation. 


13. y" = y' 14. Y= yy 3) 


In Problems 15 and 16 interpret each statement as a differen- 
tial equation. 


15. On the graph of y = (x) the slope of the tangent line at 
a point P(x, y) is the square of the distance from P(x, y) to 
the origin. 


16. On the graph of y = f(x) the rate at which the slope 
changes with respect to x at a point P(x, y) is the nega- 
tive of the slope of the tangent line at P(x, y). 


17. (a) Give the domain of the function y = ST. 


(b) Give the largest interval J of definition over which 


y =x? is solution of the differential equation 
3xy’ — 2y = 0. 
18. (a) Verify that the one-parameter family y? — 2y = 


x? — x + cis an implicit solution of the differential 


equation (2y — 2)y’ = 2x — 1. 
(b) Find a member of the one-parameter family in 
part (a) that satisfies the initial condition y(0) = 1. 


(c) Use your result in part (b) to find an explicit 
function y = (x) that satisfies y(0) = 1. Give the 
domain of the function dh. Is y = d(x) a solution of 
the initial-value problem? If so, give its interval J of 
definition; if not, explain. 


19. Given that y = x — 2/x is a solution of the DE xy’ + 
y = 2x. Find xo and the largest interval J for which y(x) is 
a solution of the first-order IVP xy’ + y = 2x, y(xo) = 1. 


20. Suppose that y(x) denotes a solution of the first-order 
IVP y’ = x? + y?, y(1) = —1 and that y(x) possesses 
at least a second derivative at x = 1. In some neigh- 
borhood of x = use the DE to determine whether 
y(x) is increasing or decreasing and whether the graph 
y(x) is concave up or concave down. 

21. A differential equation may possess more than one fam- 
ily of solutions. 

(a) Plot different members of the families 
vs dt zs +c and y= (x) = =x + c. 

(b) Verify that y= db) and y= (x) are two 
solutions of the nonlinear first-order differential 
equation (y')? = 4x?. 

(c) Construct a piecewise-defined function that is a 
solution of the nonlinear DE in part (b) but is not a 
member of either family of solutions in part (a). 


22. What is the slope of the tangent line to the graph of a 
solution of y’ = 6V/y + 5x that passes through (—1, 4)? 


In Problems 23-26 verify that the indicated function is a 
particular solution of the given differential equation. Give an 
interval of definition 7 for each solution. 

23. y” + y =2cosx-— 2sinx; y=xsinx + xcosx 


24. y” + y= secx; y= xsinx + (cos x)In(cos x) 
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25. x?y"+ xy’ +y=0; y= sin(nx) 
26. x*y" + xy’ + y = sec(In x); 
y = cos(In x) In(cos(In x)) + (ln x) sin(ln x) 
In Problems 27-30, y = cje™ + coe-* — 2x is a two- 
parameter family of the second-order DE y” — 2y’ — 3y = 


6x + 4. Find a solution of the second-order IVP consisting 
of this differential equation and the given initial conditions. 


27. y (0) = 0, y'(0) =0 28. y (0) = 1, y'(0) = —3 
29. yD=4y'(D=—-2 30. y(-l)=0,y'(-D =1 


31. The graph of a solution of a second-order initial-value 


problem d*y/dx? = f(x, y, y’), Y2) = yo y'(2) = yn is 
given in Figure 1.R.1. Use the graph to estimate the val- 
ues of yo and vu, 


a5 +: 


FIGURE 1.R.1 Graph for Problem 31 


33. 


34. 


CHAPTER 1 IN REVIEW e 33 


. A tank in the form of a right-circular cylinder of radius 
2 feet and height 10 feet is standing on end. If the tank 
is initially full of water and water leaks from a circular 
hole of radius H inch at its bottom, determine a differen- 
tial equation for the height h of the water at time t. 
Ignore friction and contraction of water at the hole. 


The number of field mice in a certain pasture is given by 
the function 200 — 10t, where time f is measured in 
years. Determine a differential equation governing a 
population of owls that feed on the mice if the rate at 
which the owl population grows is proportional to the 
difference between the number of owls at time f and 
number of field mice at time t. 


Suppose that dA/dt = —0.0004332 A0) represents 
a mathematical model for the radioactive decay of 
radium-226, where A(t) is the amount of radium (mea- 
sured in grams) remaining at time ¢ (measured in years). 
How much of the radium sample remains at the time t 
when the sample is decaying at a rate of 0.002 gram per 
year? 
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ANSWERS FOR SELECTED 


ODD-NUMBERED PROBLEMS 


EXERCISES 1.1 (PAGE 10) 


1. linear, second order 3. linear, fourth order 
5. nonlinear, second order 7. linear, third order 
9. linear in x but nonlinear in y 
15. domain of function is [—2, œ); largest interval of 
definition for solution is (—2, ©) 
17. domain of function is the set of real numbers except 
x = 2 and x = —2; largest intervals of definition for 
solution are (~œ, —2), (—2, 2), or (2, ©) 


t=] 
19. X = a ` defined on (—%, In 2) or on (In 2, œ) 
SE 
27. m= —2 29. m= 2,m=3 31. m = 0,m = —1 
33. y=2 35. no constant solutions 


EXERCISES 1.2 (PAGE 17) 
1. y = 1/0 — Aen 


3. y= 1/@? — 1); (1, %) 

5. y = 1/(x* + 1); (~, %) 

7. x= —cost+ 8 sin t 

9. x=Bcost+isint 11. y=3e'—se% 
13. y = äer! 15. y=0,y=x3 


17. half-planes defined by either y > 0 or y < 0 
19. half-planes defined by either x > 0 or x < 0 
21. the regions defined by y > 2, y < —2,or —2 <y < 2 
23. any region not containing (0, 0) 
25. yes 
27. no 
29. (a) y=cx 
(b) any rectangular region not touching the y-axis 
(c) No, the function is not differentiable at x = 0. 
31. (b) y = 1/(1 — x) on (—%, 1); 
y = —1/@ + 1) on(-1, %); 
(c) y = 0 on (—%, %) 


EXERCISES 1.3 (PAGE 27) 


dP dP 
1. — = kP + r,— = kP -r 
dt dt 
dP o 
3. a =k,P— kb 
7 Oy (1000 — x) 
Va x x 
9 EE 
` dr 100 ` 
dA 7 dh CT 
11. + A= 13. — = — — Vh 
dt 600-t $ S dr 450 V" 


di : dv 

15. L— + Ri = E(t) 17. m— = mg — kv’ 
dt dt 
dx dr gR _ 

19. mpz = —kx T CR 


dA dx 
23. — = k(M — A), k > 25. — + kx = r,k > 
3 Es ( ) 0 ES x=r 0 


dy -x+ Ve +y 


27. 
dx y 


CHAPTER 1 IN REVIEW (PAGE 32) 


1. 2 = 10y 3. y" + Ry =0 
dx 

5, y"—2y' ty =0 7. (a), (d) 

9. (b) 11. (b) 


13. y = cı and y = c2e’, cı and c2 constants 

15. y' =x? + y? 

17. (a) The domain is the set of all real numbers. 
(b) either (—~, 0) or (0, œ) 

19. For xọ = —1 the interval is (—~, 0), and for xọ = 2 the 
interval is (0, ©). 


=x, x<0 


25. (0, ~) 
29. y= Zei"? oe 8th = 2x, 


31. yo = —3,y1 =0 


dP 
33. ES = k(P — 200 + 10ż) 


This page contains answers for this chapter only. 


ANS-1 


Copyright 2009 Cengage Learning, Inc. All Rights Reserved. May not be copied, scanned, or duplicated, in whole or in part. 


